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1. Introduction 

The purpose of this paper is to discuss the consistency, methodology and accuracy of 
a very long calculation in perturbative quantum gravity. Although such an exposition is 
necessary if the result is to win acceptance, we should best begin by explaining the mo- 
tivation. This labor was undertaken to check the suggestion [1,2] that corrections from 
the infrared of quantum gravity may extinguish inflation without the need for a special 
inflaton field. It will be seen that this is also a proposal for solving the problem of the 
cosmological constant. In our scheme the cosmological constant is not unnaturally small, 
it only appears so today on account of screening by infrared effects in quantum gravity. 
Inflation begins in Hubble-sized patchs of the early universe over which the local temper- 
ature has fallen enough for the cosmological constant to dominate the stress-energy. After 
a few e-foldings the temperature of such a patch is sufficiently low to permit long range 
correlations, and the infrared screening effect begins to build up. This build-up requires 
a long time because gravitational interactions are naturally weak; they can only become 
significant through causal and coherent superposition over the past lightcone. Effective 
screening is therefore delayed until an enormous invariant volume has developed within 
the past lightcone of the observation point to the onset of inflation. This is why inflation 
lasts long enough to explain the homogeneity and isotropy of the observed universe. 

Note that no fine tuning is necessary for our scheme, nor do we require new matter 
fields or even new gravitational interactions. Of the phenomenologically viable quanta, 
our mechanism is unique to gravitons. The other known particles are either massive - 
which precludes coherent superposition — or else they possess conformal invariance on the 
classical level — which means they cannot exploit the enormous invariant volume in the 
past lightcone of a conformally flat, inflating universe. 

Note also that causality, and the physically motivated initial condition of coherent 
inflation over a finite spatial region, preclude sensitivity to global issues such as whether 
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or not the full de Sitter manifold is used. For convenience we worked on the manifold 
T 3 x 3?. What we computed is the expectation value of the invariant element in the 
presence of a homogeneous and isotropic state which is initially free de Sitter vacuum: 





g^(t, x) dx^dx v Oj = -dt z + a z (t) dx ■ dx . (1.1) 

Inflation redshifts the temperature to zero so rapidly that we simply used zero temperature 
quantum field theory. We were able to use the Lagrangian of conventional general relativity: 

C = -J-(R-2A)V^ , (1-2) 



16irG 

because our infrared mechanism is insensitive to the still unknown ultraviolet sector of 
quantum gravity. We assumed only that the scale of inflation M ~ (A/G) 1 / 4 is at least a 
few orders of magnitude below the Planck mass Mp\ ~ G^ 1 / 2 ); 

M < 1(T 3 M P1 . (1.3) 

In this case ultraviolet modes should have plenty of time to reach their natural equilibria 
as they redshift down to scales at which the dynamics is described by quantum general 
relativity. Our mechanism can be shown to derive from modes whose physical wavelength 
is about the Hubble radius at the time that they contribute most strongly. It is also worth 
pointing out that our mechanism is an inherently quantum mechanical effect, deriving 
ultimately from the gravitational interaction between the zero point motions of the various 
modes. It in no way conflicts with the classical and semi-classical stability of locally de 
Sitter backgrounds [3]. 

We inferred the physical rate of expansion from the effective Hubble parameter: 

H eS (t) ee j f ln(a) . (1.4) 
The first secular effect occurs at two loops and has the form: 

H eS (t) = H^l- I (Htf + (subdominant) + 0(k 6 )| , (1.5) 



where k = 16nG is the usual loop counting parameter of quantum gravity and 3H 2 = A 
is the bare Hubble constant. We were also able to show that the dominant contribution 
at £ loops has the form: 

-# (kH) 2£ (Ht) £ . (1.6) 

This permits the following estimate of the number of e-foldings needed for the extinction 
of inflation: 

N = Ht ~ (kH)- 2 > 10 12 . (1.7) 

Of course perturbation theory is no longer reliable when quantum corrections become of 
order one so the valid conclusion is that quantum gravitational effects slow the rate of 
inflation by an amount which becomes non-perturbatively large at late times. 

This completes our discussion of why we undertook the project and what we got. A 
more complete discussion of the physical consequences can be found elsewhere [4]. The 
remainder of this paper is devoted to an explanation of how we obtained the result and 
why we believe it is accurate. Section 2 deals with basic issues of the formalism. Crucial 
points here are the distinction between "in" -"out" matrix elements and expectation values, 
the physical significance of our ultraviolet regularization, the demonstration that H e g(t) 
is independent of the choice of gauge, the reason why infrared logarithms are almost 
inevitable, and a derivation of our bound (1.6) for the dominant corrections at I loops. 
Sections 3 and 4 deal with our specific process. Section 3 is devoted to the class of diagrams 
which involve only 3-point vertices, while Section 4 describes the dominant diagram which 
has a 4-point vertex and a 3-point vertex. Section 5 assembles the final result and discusses 
the many accuracy checks. We close the section, and the paper, with a brief consideration 
of the issues pertaining to the case of negative A. 
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2. Basic Issues 

2. 1 The Apparatus of Perturbation Theory. 

We take the onset of inflation to be t = 0, and we work perturbatively around the 
classical background: 

a c lass(*) = exp(tft) . (2.1) 

It is simplest to perform the calculation in conformally flat coordinates, for which the 
invariant element of the background is: 

-dt 2 + al lass (t) dx ■ dx = Q 2 (-du 2 + dx ■ dx) , (2.2a) 

O = {Hu)~ l = exp(#£) . (2.26) 

Note the temporal inversion and the fact that the onset of inflation at t = corresponds 
to« = H -1 . Since t — > oo corresponds to u — > + , and since the spatial coordinates of T 3 
fall within the region, — \H~ l < x l < \H~ X , the range of conformal coordinates is rather 
small. This is why a conformally invariant field — whose dynamics are locally the same 
as in flat space, except for ultraviolet regular ization — cannot induce a big infrared effect. 

Perturbation theory is organized most conveniently in terms of a "pseudo-graviton" 
field, ipfj,u, obtained by conformally re-scaling the metric: 

guv = 9uv = ^ (j) [W + Kip^ v ) . (2.3) 

As usual, pseudo-graviton indices are raised and lowered with the Lorentz metric, and 
the loop counting parameter is k 2 = 167rG. After some judicious partial integrations the 
invariant part of the bare Lagrangian takes the following form [5]: 

Anv = \f~9 Z 9 af3 9 Pa 9^ ^ 2 h^ap,u ^ua,/3 ~ \^a0,p + i^a/3 lP ^uv,a 

- \^ap,a Tpfai] -\^9 9 pa 9^ (0% <W V„ a • (2.4) 
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Gauge fixing is accomplished through the addition of — \rf v F^Fy, where [5]: 

(V« ,« - \^ + 2 ^ (lnft), «) • 
The associated ghost Lagrangian is [5]: 



(2.5) 



£ghost = -n 2 V*" 



g P H d v + g pv <9 M + 3 W + 2g liV (In ft) 



+ (ft 2 cl^) 77^ fe p a CT + + g pa (In ft) J . (2.6) 



In our gauge the pseudo-graviton kinetic operator has the form: 



D pa 



W S ^ - \^ ^ - D A - t {fl S u) [p H>) D B + t^ v tH° D C , (2.7) 



where parenthesized indices are symmetrized. Two notational conventions reflect the fact 
that the direction is special. First, we define as: 



(2.8a) 



(Recall from (1.1) that our metric is spacelike.) Second, we define barred tensors to have 
their natural zero components nulled, for example: 



(2.86) 



Note also that = ft(<9 2 + Jr)ft is the kinetic operator for a massless, minimally coupled 
scalar and = = ft <9 2 ft is the kinetic operator for a conformally coupled scalar. 

The zeroth order action results in the following free field expansion [6] : 
ip^(u,x) = (J^J +H 3 Yl \^fiv{u,x;k,Xj a(k, A) + (u, x; k, Xj a\k,X) 



The spatial polarizations consist of "A" modes: 

Hu 



^Ijlv (u, x; fc, A j = —^L (l + -j^j exp ik(u - -gj + ik ■ x e^(k, A) 



(2.9) 

VA e A 
(2.10a) 



while the space-time and purely temporal polarizations are associated, respectively, with 
"B" and "C" modes: 



jj,v \ u i x \ k, A ) — 



Hu 



exp 



ik (u- i) +ik-x e^(k,\) ; VA e S, C (2.106) 



In LSZ reduction one would integrate against and contract into ^^(-u, x; k, A) to insert 

— * 

and "in" -coming graviton of momentum k and polarization A; the conjugate would be 
used to extract an "out" -going graviton with the same quantum numbers. The zero modes 
evolve as free particles with time dependences 1 and for the A modes, and u and v? for 
the B and C modes. Since causality decouples the zero modes shortly after the onset of 
inflation, they play no role in screening and we shall not trouble with them further. 

We define |0) as the Heisenberg state annihilated by a(k 7 A) — and the analogous ghost 
operators — at the onset of inflation. We can use this condition and expansion (2.9) to 
express the free pseudo-graviton propagator as a mode sum [7] : 



^Apcr] (x;x) = ^o|t|^(x) ip pa (x') 



free: 



(2.11a) 
(2.116) 



Note that the convergence factor e e l I fc l I serves as an ultraviolet mode cutoff. Although 
the resulting regularization is very convenient for this calculation, its failure to respect 
general coordinate invariance necessitates the use of non-invariant counterterms. These 
are analogous to the photon mass which must be added to QED when using a momentum 
cutoff. Just as in QED, these non-invariant counterterms do not affect long distance 
phenomena. 

Because the propagator is only needed for small conformal coordinate separations, 
Ax = \\x ' — x\\ and Au = u' — u, the sum over momenta is well approximated as an 
integral whose lower limit is the momentum k = H of the longest wavelength. When this 
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is done the pseudo-graviton propagator becomes [6] : 
d 3 k (H 2 uu' 



exp 



H 2 (l + ik\Au\) 



+ 



2k?> 



— ik\ Au\ + ik ■ (x ' — x) — ek 
exp 



2* M (p 8 v a) - r,^ v pa 



ik\ Au\ + ik ■ (x ' - x) - ek] \28^ p 8 u a) - 2rJ^ 



(2.12a) 



so that: 

A^] (a; x') 



[IV 



2u'u 



87T 2 \ Ax 2 - Au 2 + 2ie\ Au\ + e 2 
-In 



28^ C } - Vp, V pa 



H 2 (Ax 2 - Au 2 + 2ie\Au\ + e 2 ) 



28 } p 5^ 



2V*T \ (2-126) 



The same approximation gives the following result for the ghost propagator: 

H 2 { 2u'u 



(x;x') 



8tt 2 I Ax 2 - Au 2 + 2ie\ Au\ + e 2 
-In 



H 2 (Ax 2 - Au 2 + 2ie\Au\ + e r 



(2-13) 



The decoupling between tensor indices and the functional dependence upon spacetime - 
and the simplicity of each — greatly facilitates calculations. It is convenient to identify as 
the "normal" and "logarithmic" propagator functions as iA_/y and iA^ respectively: 

H 2 2uu' 



iAn(x, x) 



8tt 2 Ax 2 - Au 2 + 2ie\Au\ + e 2 



H 2 ( Au 2 - Au 2 + 2ie\Au\ + e 2 



H 2 

iAr(x,x) = — o In 
LK ' 8tt 2 



In this notation we can write the pseudo-graviton and ghost propagators as: 



M „A^] (x; x') =iA N (x; x') [28„ {p 8j ] - ^ V P ° 

-iA L (x;x') \28} p 8j ) -2rj ul/ rj p( > 



(2.14a) 
(2.146) 



iiA v 



(x; x) = iA N (x, x) r]^ - iA L (x; x) 7]^ 



(2.15a) 



(2.156) 
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2.2 "In" -"Out" Matrix Elements and the S- Matrix. 

Perturbative quantum field theory is usually formulated to give "in" -"out" amplitudes 
and S-matrix elements. These quantities are not well suited for our study because they re- 
quire specification of the vacuum at asymptotically early and late times, which is precisely 
what we wish to determine. However, they can at least be used to negate the hypothesis 
that the vacuum of an initially inflating universe suffers only perturbatively small correc- 
tions. The procedure is simply to assume the "in" and "out" vacua are both free de Sitter, 
and then do the computation. If the hypothesis is correct, the result should be free of 
infrared divergences. 

What one actually finds is that both "in"- "out" amplitudes [1,2] and S-matrix elements 
[7] are infrared divergent. Since even 3-particle tree amplitudes are affected, there is no 
possibility for solving the problem by summing degenerate ensembles. There is simply 
nothing of lower order in perturbation theory that could cancel the problem. 

Why this happens can be readily understood from the previous section. S-matrix 
elements consist of wavefunctions integrated against interaction vertices, which are linked 
by propagators. From (2.10a) we see that physical wavefunctions become constant at late 
times; while (2.14-15) shows that non-coincident propagators remain of order one as u — > 0. 
But (2.4) reveals that vertices blow up at late times. A more physical way of understanding 
the phenomenon is that although the coordinate momentum of a graviton is unchanged by 
time evolution, its physical momentum is redshifted to zero. Since all gravitons of fixed 
coordinate momentum approach the same physical momentum, their interaction becomes 
infinitely strong at late times. 

The correct interpretation of these infrared divergences is that the "in" vacuum is 
infinitely far from the "out" vacuum. Stated differently, the vacuum of an initially inflating 
universe suffers non-perturbatively large corrections at late times. However, it does not 
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follow that these corrections must depend upon time in the same way that the "in" -"out" 
divergences depend upon some arbitrarily chosen infrared cutoff. To obtain a quantitative 
result one must actually follow the evolution. 

2.3 How to Compute Expectation Values. 

The perturbative rules for calculating expectation values of operators in field theory 
were developed by Schwinger [8] and have been adapted to our particular problem and 
initial conditions in [2]. They are quite similar to the usual rules for computing "in" -"out" 
matrix elements since the propagators and vertices utilized are simple variations of the 
usual ones. The main idea is to evolve forward from the initial state with the action 
functionally integrated over the dummy field ipffi and, then, to evolve back to the initial 
state using the conjugate action functionally integrated over the dummy field ipjn) . 

The Feynman rules are simple. An external line may be chosen as either "+" or "— " 
but one does not sum over both possibilities. Vertices have either all "+" or all "— " lines, 
and one sums the two possibilities for each vertex. The "+" vertices are identical to those 
of the "in"- "out" formalism, while the "— " vertices are complex conjugated. Propagators 
can link fields of either sign. All four possibilities can be obtained from (2.12a-b) by 
effecting the following substitutions for the quantity (Ax 2 — Au + 2ie\Au\ + e 2 ) in (2.14): 



(+ +) = 


=> (Ax- 


\Au\ + ie) (Ax + \Au\ - ie) , 


(2.16a) 


(- +) = 


=>> (Ax - 


Au + ie) (Ax + Au- ie) , 


(2.166) 


(+ -) = 


(Ax- 


Au - ie) (Ax + Au + ie) , 


(2.16c) 


(- -) = 


^ (Ax - 


|Au| - ie) (Ax + \Au\ + ie) . 


(2.16d) 



Two properties of Schwinger's formalism have significance for our computation. First, 
the expectation value of a Hermitian operator such as the metric must be real. Second, 
the interference between "+" and "— " vertices results in complete cancellation whenever 
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an interaction strays outside the past lightcone of the observation point. This means that 
no infrared regularization is needed. The infrared divergences of the "in"- "out" formalism 
come from interactions in the infinite future, and these drop out of "in" -"in" calculations. 
In their place one expects growth in the observation time. 

2.4 Ultraviolet Regularization. 

The simplest way of regulating the ultraviolet is by keeping the parameter e non- 
zero in the "++" propagators (2.14-15) and their variations (2.16). Consideration of the 
mode sums (2.11) reveals e^ 1 to be an exponential cutoff on the coordinate 3-momentum. 
This would not be a very natural technique for an "in"- "out" calculation because there 
is nothing unique about the t = surface of simultaneity. However, this surface has a 
special significance in our modified "in" -"in" computation: it is where the initial state 
is defined, and it is the point at which interactions begin. In the context of expectation 
values, our method corresponds to weighting the usual Fock space inner product with a 
factor of exp(— ^e|| k\\) for each creation and annihilation operator. In other words, time 
evolution is unaffected by our method, only the inner product on the initial value surface 
changes. 

The point just made is crucial because spurious time dependence can be made to 
reside on the ultraviolet regularization parameter. Consider, for example, a logarithmic 
divergence of the form ln(e). Had we instead regulated by replacing e everywhere with 
eHu, then we might have claimed to see a secular infrared effect: 

In(eHu) = ln(e) - Ht (2.17) 

which is in fact of purely ultraviolet origin. We emphasize that the possibility for this 
sort of delusion has nothing to do with the perturbative non-renormalizability of quantum 
general relativity. One can see it even in the flat space, massless <f> model which was our 
original paradigm for relaxation [2] . Nor does the phenomenon signify any fuzziness in the 
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distinction between infrared and ultraviolet. The unambiguous signal for an ultraviolet 
effect is propagators at or near coincidence. 

The time dependence of the inflating background has led to much confusion about 
ultraviolet regularization. Some researchers have thought it more natural to employ a 
mode cutoff which is invariant with respect to the background geometry. (Note that this 
method is no more invariant than ours with respect to the full geometry.) This amounts 



to replacing our parameter e with eV H 2 uu' . One consequence is that the coincidence 
limit of an undifferentiated propagator which has been so regulated grows linearly in the 
co-moving time. By first taking the coincidence limit of our propagator and then replacing 
e with eHu: 

H 2 { 2u'u 



[IV 



A pa 



(x; x) 



8n 2 



-In 



H 2 e 2 



8tt 2 1 H 2 e 2 



25^ 6 v a) 



fW 

pa 



(2.18o) 



+ (2Ht - In 



H 2 e 2 



2S} P 5° } - 2r) uy rP 



(2.186) 



we see that this growth is identical to that displayed in (2.17), and hence of purely ultra- 
violet provenance. One gets the same result with background-invariant point splitting. 

An erroneous argument is sometimes given that this spurious time dependence from 
the ultraviolet is actually a reliable effect from the infrared. One first notes that on 
the propagator would be our integral expression (2.12a), but without the lower bound 
at k = H. The second term of the integrand grows rapidly enough near k = to give 
an infrared divergence. According to the fallacious argument, one should regulate this 
infrared divergence by cutting off the integration at k = (uu')~' ! 2 . Of course this replaces 
the factor of H 2 inside the logarithm term of our propagator with (uu f )^^, and one can 
see from (2.18a) that the coincidence limit then exhibits the same linear growth as (2.18b). 
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We stress that the introduction of a time dependent infrared cutoff prevents the prop- 
agator from even inverting the kinetic operator. The rationale behind the cutoff is that 
coordinate momenta smaller than k = correspond to physical wavelengths which have 
redshifted beyond the Hubble radius and should therefore decouple. This is correct physics 
but faulty mathematics. Modes indeed decouple when their physical wavelengths redshift 
beyond the Hubble radius, but this is accomplished by the causality of interactions in 
Schwinger's formalism, not by the ad hoc imposition of an infrared cutoff on the naive 
mode sum. The mode sum is actually regulated by the physically motivated restriction of 
coherent de Sitter vacuum to a spatial patch of finite extent. For us this was the Hubble ra- 
dius, but any value would serve. The key point is that the infrared cutoff on the free mode 
sum is not time dependent. The only time dependence which appears in the coincidence 
limit of the propagator got there through a poor choice of the ultraviolet regulator. 

Motivated by the faulty argument — for which we emphasize that he bore no respon- 
sibility — Ford [9] proposed a very interesting relaxation mechanism of which ours is, in 
some ways, a mirror image.* Among other differences, the most important diagrams for 
Ford were those with a coincident propagator attached to two legs of the same vertex. 
These are precisely the least important ones for us. 

2.5 The Threshold for a Late Time Effect. 

We found it convenient to compute the amputated expectation value of the pseudo- 
graviton field and then attach the external line by solving an ordinary differential equation. 
The homogeneity and isotropy of the dynamics, and of our initial state, allow us to express 
the amputated 1-point function in terms of two functions of u: 



The full 1-point function must have the same form, although with different coefficient 



We thank M. B. Einhorn for bringing Ford's work to our attention. 




(2.19) 
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(2.20) 



functions: 

(o| k% v {u, x) |o) = A(u) fj^ + C{u) S°S V ° . 

Comparing these two expressions and taking account of the kinetic operator (2.7), we see 
that A(u) and C(u) can be expressed in terms of a(u) and c(u) using the retarded Green's 
functions for the massless, minimally coupled and conformally coupled scalars: 



A(u) = -4 Of a (u) + G r S 



net 



3a + c 



(«) , 



ret 



C(«) = 



3a + c 



(«) 



(2.21a) 



(2.216) 



Although it is simple to give integral expressions for these Greens functions, the more 
revealing form is to work them out for arbitrary powers of u [2]: 



^ret 



u~ A (Hu) C (u) 
u- 4 (Hu) C ] (u) 



sc(i - K) 
if 2 



{Hup - (i - k) - k(#«) J 



(2.22a) 



-(£Tu)^+2(l-iC)^«-(l-C)(^w) 2 [• ( 2 -226) 



2(1 - o(i -ho 

Note that ( = constitutes a sort of threshold. For larger values of ( the late time limits 
of A(u) and C{u) approach constants, whereas smaller values of ( give functions which 
grow as u approaches zero. 

To obtain H e ^{t) we compare with the invariant element in co- moving coordinates 

-dt 2 + a?(t) dx-dx = tt 2 {-[l-C(u) du 2 + 1 + A(u)j dx ■ dx } . (2.23) 

Substituting into the definition (1.4) gives a result which is true even beyond perturbation 
theory: 



H eS (t) = 



H 



! _ 2" luMu) 



(2.24) 



v/l - C(u) I l + A(u) 
With relations (2.21) and (2.22), this shows that a(u) and c(u) must grow faster than -u~ 4 
if quantum corrections are to overwhelm the classical result of H c g(t) = H at late times. 
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2.6 Gauge Independence of H e g(t). 

In demonstrating the gauge independence of H c ^(t) it is useful to recall that any 
quantity can be made invariant by defining it in a unique coordinate system. In this case 
we can also exploit the homogeneity and isotropy of our special state, although we feel 
confident that a satisfactory definition of the effective Hubble constant could be found for 
more general initial states. Since H c g(t) was defined in co-moving coordinates, for a state 
whose scale factor obeys the initial conditions: 

a(0) = l , ^ = ff; (2.25) 

our first step in any new coordinate system would be to restore these conditions by per- 
forming an appropriate coordinate transformation. We can therefore limit ourselves to 
consideration of gauge changes — > F fl + 5F^ which preserve homogeneity, isotropy, and 
the initial conditions. Any such change can be parametrized as follows: 

SFfj, = 5® (p(u) , (2.26) 

where preserving the initial condition requires (p(H~^) = 0. 

Changing our gauge condition (2.5) by (2.26) induces the following 1-point interaction: 

6 ( K ~i J d A x F^F^j = J d A x F ft - Vu) (2.27a) 

= J d 4 x j^oo <f,0 + \i> <P,0 + I ^00 V?} • (2-276) 

The coefficient functions of the amputated 1-point function therefore acquire the following 
variations: 

8a(u) = | <p(u) , 5c(u) = | <p(u) + ^ <p(u) . (2.28) 

Now consider variation by a general power, minus a constant to enforce the initial condition: 

(p(u) = (Hu) C - H 3 . (2.29) 
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The induced variations on a(u) and c(u) are: 



5a{u) = - (C - 3) u- 4 (Hu) C , (2.30a) 

5c(u) = £ (C + 1) (#w) C - =^==- • (2.306) 

The coefficient functions of the unamputated 1-point function suffer the following varia- 
tions: 

^)— 2 {-f^ + ^-(^)f + ^}, (,3 l0 ) 

, cw = 2 ^{_0_(£^)^ + (^)!}. ( , 316) 

Prom (2.24) we see that the variation of H e ^{t) is: 

MfefKf) = H |~ « ^ M(«) + i , (2.32) 

and substitution of (2.31) gives 5H e ff(t) = 0. The proof is completed by noting that 
a general variation of the required form can be built by superposing terms of the form 
(2.29). 

2. 7 The Genesis of Infrared Logarithms. 

It turns out that perturbative corrections to the amputated 1-point function can grow 
no faster than -u~ 4 times powers of \n(Hu) [2,10]. Since we will need to extend the ar- 
gument, it is useful to begin by restating it. Recall first that k has the dimensions of 
length, while H goes like an inverse length, and the amputated 1-point function has the 
dimensions of an inverse length squared. Diagrams which contribute to the amputated 
1-point function at t loops can involve up to 2£ — 1 3-point interactions and 3£ — 2 propa- 
gators. Each 3-point vertex contributes a constant factor of kH~ 2 , while each propagator 
contributes a factor of H 2 . Including the single factor of k in our definition (2.19), we see 

Of P — 

that the £-loop contribution consists of k H times a function of H, u and e which has 
the dimensionality (length) ~ 4 . 
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Since the ultraviolet regularization parameter goes to zero after renormalization, we can 
ignore positive powers of e. The constant H can enter through undifferentiated iA^, and 
through the limits of integration on the conformal time integrals. Neither can contribute 
negative powers of H. This is obvious for iA^. To see that factors of H from the limits 
cannot contribute negative powers, note that performing the spatial integrations leaves 
2£ — 2 conformal time integrations of an integrand which falls off like g -2 ^ 6 if a ll the 
conformal times are scaled by the common factor S. Since the integrand can contain at 
most two isolated factors of the conformal time integrations must converge even if the 
factors of H~^~ in their upper limits are taken to infinity. Hence there can be no negative 
powers of H. 

To complete the argument consider a dimension (length)" 4 function of H, u and e which 
contains no negative powers of H or positive powers of e. We can write any such function 
as -u™ 4 times a combination of the dimensionless parameters Hu and ew _1 . Since there can 
be no negative powers of the first parameter, or positive powers of the second, the strongest 
growth possible as u — > + is w~ 4 times powers of ln(Hu) and ln(ett _1 ) = ln(He) — ln(Hu). 

Because the amputated 1-point function must grow faster than w~ 4 in order for there 
to be a significant effect at late times, the appearance of infrared logarithms is essential for 
our scheme of relaxation. The preceding argument can be extended to show that infrared 
logarithms are just about inevitable if one accepts that there are logarithmic ultraviolet 
divergences. Since e is a dimensionful quantity, logarithmic ultraviolet divergences can 
only come in the form of ln(ett _1 ) and ln(He). As before, there are just two sources of 
H dependence: an undifferentiated iA^ and the upper limits H of the conformal time 
integrations. The latter cannot provide factors of ln(H) for the same reason it gives no 
negative powers of H: the conformal time integrands fall off rapidly enough for large 
conformal times to make them converge even when the upper limits become infinite. The 
logarithm part of an undifferentiated propagator can indeed provide a factor of In (if), but 
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there are never enough such terms to pair with all the factors of ln(e). At £ loops one has 
to expect £ factors of ln(e) from ultraviolet divergences, but there is a maximum of £ — 1 
logarithms from undifferentiated iA^'s. * Hence at least one of the ultraviolet logarithms 
must come in the form ln(eit — *). 

We emphasize that the association we have exploited between infrared logarithms and 
logarithmic ultraviolet divergences in no way implies that factors of In (it) are of ultraviolet 
origin. Some of them originate as factors of ln(Hu) from the logarithm parts of undiffer- 
entiated propagators. ** The physical origin of such terms is the increasing correlation 
of the free graviton vacuum as inflation proceeds. This effect is the casual analog of the 
problem discussed in sub-section 2.4 when one assumes the existence of coherent de Sitter 
vacuum over an infinite surface of simultaneity. 

Even the factors of ln(*u) which originate as h^e-u^ 1 ) are infrared effects. They derive 
from the coherent superposition of interactions throughout the invariant volume of the 
past lightcone. In the absence of a mass, ultraviolet divergences must always be associated 
with the infrared in this way. See, for example, the result we obtained for massless, 4? 
theory in flat space [2]. 

2.8 A Bound on the Number of Infrared Logarithms. 

We argued in the preceding sub-section that infrared logarithms are all but inevitable, 
and that they derive jointly from the infrared regions of loop integrals which harbor log- 
arithmic ultraviolet divergences and from an undifferentiated iA^. However, we did not 
explain how many of the £ — 1 undifferentiated propagator logarithms can go to reinforce 
the i ultraviolet logarithms that one expects in an £-loop diagram. The answer turns out 

See the end of sub-section 2.9 for a proof. 
The reader should not be misled by our argument that there are fewer intrinsic factors of hn(Hu) than of 
ln(eu _1 ). We noted this only to rule out the remote possibility that the two factors come in pairs so as to 
cancel the infrared logarithms: ln(eu _1 ) + \n(Hu) = ln(He). In fact the two sources of infrared logarithms 
tend to add. 
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to be "none:" there can be at most I infrared logarithms at I loops. All the undifferenti- 
ated factors of iA^ do is to sometimes make up the difference for diagrams which would 
otherwise contribute less than £ infrared logarithms. 

To see why, consider the way a single undifferentiated propagator might enhance the 
number of infrared logarithms. There can be I logarithmic ultraviolet divergences at I 
loops so what we are looking for is: 

ln*(eu _1 ) x \n(Hu) = ln^ew" 1 ) x ln(H) + ln*(eu _1 ) x ln(u) . (2.33) 

In other words, if we replace the propagator in question with ln(if), the resulting diagram 
must contribute £ logarithmic ultraviolet divergences. But replacing a propagator by a 
constant cuts the leg and results in a diagram containing only I — 1 loops. Such a diagram 
can contribute at most £ — 1 factors of ln(ew _1 ). Therefore, an undifferentiated iA^ cannot 
increase the number of infrared logarithms in an £-\oop diagram beyond £. 

It is worth mentioning that this argument applies to diagrams, not to portions of 
diagrams. However, the only practical way of performing the computation was to break 
each diagram up into many pieces. At two loops some of these pieces produce three 
infrared logarithms; it is only their sums which are limited to two. On the first run of this 
calculation we did not appreciate that the triple log terms must cancel. This argument was 
only discovered after noting what was to us then, a surprising and disturbing cancellation. 
In retrospect the cancellation stands as a powerful check on the consistency and accuracy 
of our work. 

2.9 Comments on One and Two Loops. 

The first quantum corrections to the amputated 1-point function are the one-loop 
graphs shown in Fig. 1. Since the external line has been amputated, there is no integration 
over the single interaction vertex at x^ 1 = (u, x). One obtains just the coincidence limit of 
a pseudo-graviton or ghost propagator, acted upon by the appropriate vertex operator. It 
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turns out that there can be no undifferentiated "logarithmic" parts iA^, but even if there 
were, they would only contribute factors of \n(H e) . One needs to integrate something over 
a large invariant volume in order to see an effect. This is one way that our mechanism 
differs from Ford's [9]. 




00 



Fig. 1: One-loop contributions to the background geometry. Gravitons reside on wavy lines and ghosts 
on segmented lines. 

The one- loop graphs of Fig. 1 contribute at most terms of order w~ 2 , which is well below 
the threshold for producing an effect at late times. 




(a) (b) (c) (d) (e) (f) 



Fig. 2: Two- loop contributions to the background geometry. Gravitons reside on wavy lines and ghosts 
on segmented lines. 

The graphs which contribute at two loops are displayed in Fig. 2. Diagram (2f) is 
another ultra-local coincidence limit, so it contributes no infrared logarithms. Diagram 
(2e) contains a single free interaction vertex, but the entire diagram is canceled by the 
counterterm needed to renormalize its coincident lower loop. Daigram (2d) also contains 
a single free interaction, and an undifferentiated iA^ allows it to contribute two infrared 
logarithms. We can this the "4-3" diagram because it consists of a 4-point and a 3- 
point vertex. What we call "3-3-3" diagrams are shown in (2a-c). They contain two free 
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interaction vertices and would give two infrared logarithms if the calculation were done in 
flat space. By the argument of the preceding sub-section we know that undifferentiated 
logarithms from the propagators do not change this. 




(a) (b) 

Fig. 3: In the 4-3 diag ram (a) and the generic 3-3-3 diagram (b), solid lines can represent gravitons or 
ghosts and the arrows indicate the action of the vertex derivatives. 

Fig. 3 is useful in understanding why two-loop contributions to the 1-point function 
have at most a single undifferentiated iA^. One can see from (2.4) and (2.6) that every 
interaction vertex contains two fields which are either differentiated or else contain a 0- 
index. From (2.15) one can see that a 0-index precludes coupling to the logarithm term in 
the propagator. We call a line which is either differentiated or forced to carry a 0-index, 
"contaminated." 

To treat the case for a general loop £, note that a graph of this order can have 2£ — 1 
3-point vertices and 3£ — 2 internal propagators. One of the vertices is the external one, 
whose external leg may be contaminated, but the other vertices must each contribute two 
contaminated legs. This makes for a total number of either 4£ — 3 or 4£ — 4 contaminating 
factors (derivatives or 0-indices) to distribute among the 3£ — 2 internal propagators. The 
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smallest number of "spoiled" propagators is therefore 2£ — 1, leaving I — 1 which can 
harbor undifferentiated logarithms. The result is the same for £-loop graphs constructed 
from higher point vertices because the number of contaminating factors falls by two for 
each internal line which is lost. 

3. The 3-3-3 Diagrams 

The three 3-3-3 diagram consists of an outer 3-point vertex joined to two freely inte- 
grated 3-point vertices (see Fig. 4). 



|i V |i V [IV 




(a) (b) (c) 



Fig. 4: The tensor structure of the 3-3-3 diagrams. Gravitons reside on solid lines and ghosts on 
segmented lines. 

The fixed vertex is taken to be at = (u, x), and we can identify the locations of the two 
free vertices as x 1 ^ 1 = (u 1 ', x ') and x"^ = {u" ', x "). One can form three differences from 
these positions: 

w** = (x - x'Y , = (x' - x"Y , z» ee (x" - xf . (3.1) 

The fixed vertex is taken to be of "+" type, and there are a total of four variations when 
each of the two free vertices is summed over "+" and "— " . These variations have no effect 
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until quite late in the reduction process because they control only the overall sign of the 
diagram and the order e terms in the four propagators. However, since the variations 
interfere destructively whenever x ,fl or x"^" is outside the rather small past lightcone of 
x^, we will make no error by extending the spatial integrations from T 3 to 9ft 3 . 

The Lorentz squares of the vectors (3.1) are understood to include terms of order e 
according to the following scheme: 

up 1 = (r — d w + ie) (r + d w — ie) , f' = x — x'] (3.2a) 

z 2 = ( r "-d z +ie) (r" + d z -ie) , f" = x"-x; (3.26) 

y 2 = (||f f + f "|| -dy + ie) (||f f + r "|| +d y -ie) . (3.2c) 

Here r' = \\r'\\ , r" = ||r "|| , and the three d's are conformal time differences which depend 
upon the four "±" variations. For example, when is "+" and x llfJl is "— " we have: 

d w = \u — u\ , d z = u — u" , d y = u — u" . (3-3) 

We will eventually show that the order u < v! < u" < H^ 1 can be enforced, in which case 
the d's have the form given in Table 1. 



V(x'),V(x") 


dw 


d z 


dy 


+ , + 


u' — u 


u"-u 


u" - u' 


+ , - 


v! — u 


u-u" 


u' - u" 


- , + 


u — v! 


u"-u 


u" - u' 




u — v! 


u-u" 


u' - u" 



Table 1: The values of the d's for the four possible variations of the vertices at x and x" , 
assuming u < u <u" < H^ 1 . 



The ghost-graviton interactions can be obtained from expression (2.6): 
4host = K ° 2 {"^ ^ a ~ V ^ - ^ ^ + V ^ ^ 

+ \^ ^ - WW U a t a + UJ» J* ty 

+ §W uJ a ojW t a + " V ^ tv + fa & u V ta tu] ■ (3.4) 
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(Recall that = tj^q and ip = rj^ip^.) The associated vertex operators y^ 11 ^ are 
defined by the relation: 

£ ghost = Wi ^a 2 a 3 Vj* 1 ^ 1 " 2 " 3 (x ; <9i, <9 2 , <9 3 ) , %= 1, 10 ; (3.5) 

and are explicitly displayed in Table 2. 



# 


Vertex Operator 


# 


Vertex Operator 


1 


_ r} a 2 (a 1 ^ . ^ 


6 




2 




7 




3 


_^a2(ai Q&3 


8 




4 


— 2 T ja 2 (ai qM £a 3 


9 


l^aifli Qaz j-u 2 


5 


v a 3 (ai gfc) Q&2 


10 


2 ai/3i t a 2 £ a 3 



Table 2: Vertex operators of ghost-pseudograviton interactions without the factor of kD. 2 . 



We can read the graviton 3-point interaction off from expression (2.4): 

4£ = « o 2 { ^a,/? - r p taw - \r p ^ v 

- ^% + \r p ^ ^% + ta V + \r p ^ 

+ y ^ - ^ V>a/3,„ ^ - y a(3 V>,/3 

- JV + 2^ W V" + |^ a/3 ^ 

" !HZ V> V< „ ^ ^ + ^ * * + J Vfc„ *T t„ } . (3.6) 

It is worth noting that all but the last three terms can be checked against the previously 
published, flat space 3-point interactions [11] by merely omitting the factor of fi 2 and 
regarding ^)^ v as the graviton. 

In deriving the associated vertex operators we must account for the indistinguishability 
of gravitons. This would ordinarily be accomplished by fully symmetrizing each interaction, 
which turns out to give 75 distinct terms. For the pure graviton diagram (4a) it is wasteful 
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to first sum over the (75) 3 = 421,875 possibilities for the three vertices and then divide 
by the symmetry factor of 4 to compensate for overcounting. The more efficient strategy 
is to symmetrize the vertices only on line #1 and then sum over the interchange of lines 
#2 and #3 for only the vertex at x"^, dispensing with the symmetry factor. One saves 
over a factor of 3 this way. 

To obtain the partially symmetrized vertices one first takes any of the terms from (3.6) 
and permutes graviton #1 over the three possibilities. As an example, consider the first 
term in (3.6). Denoting graviton #1 by a breve, we obtain the following three terms: 



\ kQ 2 $ ^ + \ ^%a,(3 V> ^ + \ V>,/3 V> • (3-7) 



One then assigns the remaining two gravitons in each term as #2 and #3 in any way. For 
example, from (3.7) we could infer the following three vertex operators: 



The 43 operators which comprise the full vertex are given in Table 3. 

The various 3-3-3 diagrams can be written in terms of the vertex operators of Tables 2- 
3 and the propagators (2.15). The diagram of Fig. 4a has gravitons on both loops and 
results in the following expression: 



(3.8a) 



(3.86) 



(3.8c) 




(3.9a) 



25 



# 


Vertex Operator 


# 


Vertex Operator 


1 




22 


l^«i(«2 r7 /3 2 )/3i 9 («3 ^3) 


2 


_ T \_^a 2 (3 2 ^3^3 q(Pi ^/3i) 


23 


1^2(^3 ^AO/fe ^("i ^1) 


3 


-^rj a ^ 3 r) ai ^ d[ a2 th) 


24 


1^3(^1 ^ijft a}" 2 cf } 


4 


l^a^aa 7/ /3 2 )/3i ^( a 3 t /3 3 ) 


25 


^!)(a 3 ^3) ^2/32 


5 




26 


l^( a2 ?7 /3 2 )(ai ^1) ^a 3 /3 3 


6 


I^a 3 (ai ^/?i)/? 3 ^( a 2 ^2) 


27 


l^( a 3 ^/3 3 )(a 2 ^ai/3i 


7 


I^(«3 r y/?3)(ai ^a: 2 /3 2 


28 


^!)(a 2 ^2) ^303 


8 




29 


l^( a2 r7 /3 2 )(a 3 ^3) ^ai/3i 


9 


I^(a 2 ^/3 2 )(a 3 qPs) ^ai/3i 


30 


1^3 ^ 3 )( ai ^l) ^2/3, 


10 


l^fr a («2 r7 /?2)(a 3 ^3) 


31 


l^ai/Ji ^"2/32 ^asA d 2 . 


11 


l^ 2 /3 2 5 ("3 r7 /3 3 )(a 1 ^1) 


32 


l^aiA ^2 ^asA ^ . ^ 


12 


l^agft ^i)(a 2 ^2) 


33 


_l r} a 1 (a 2 ^/3 2 )/3i ^a 3 /3 3 ^ . ^ 


13 


-dj" 1 7^1) ( a 2 7^2) ( a 3 


34 


_l^a 2 (a3 ^/33)/3 2 ^ai/3i ^ . ^ 


14 


_ 9 («2 r7 /3 2 )(a 3 r? /3 3 )(ai Q^) 


35 


_1 ^a 3 (ai ^/3i)/3 3 ^a 2 /3 2 ^ . ^ 


15 


_^(«3 7/J 3 3 )(ai ^i)(a 2 ^2) 


36 


-Jc^ 1 r/ a 2/3 2 ^3 


16 


_l^( a 2 ^)(ai ^^(as fl#0 


37 


_l^( a 2 ^2) ^a 3 /3 3 ^ai/3i 


17 


-^df* 3 r/^)(a2 ?? /3 2 )(ai 


38 


_l^ai/3i ^a 2 (a3 ^/3 3 )/3 2 ^ . ^ 


18 


r^OK r/ /3 3 )(a 2 ^2) 


39 


_l^a 2 /32 ^"3(^1 ^/3i)/3 3 ^ . ^ 


19 




40 


l^"i)("2 ^2) (as ^3)(/?i 5 2 . dg 


20 


_l^a 2 /3 2 ^a 3 /3 3 ^/3i) 


41 


^ai)(Q!2 r? /3 2 )(a3 ^3) (A ^3 . ^ 


21 


_ 1^3/33 ^ai^i 9 («2 ^2) 


42 


a/3i) a 2 (a 3 „/3 3 )/32 
12 3 ' ' 






43 


^af 2 af 2) r/ a 3(«i ^0/33 



Xable 3: Partially symmetrized cubic pseudo-graviton vertex operators with vertex line #1 
distinguished and without the factor of kD. 2 . 
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The diagram of Fig. 4b has a ghost inner loop and gives: 
rH -i . . H -i . 43 10 

-« / du 1 / dV / du" / d\" J2 E 

Vp^ix'; d[, z[ a2 A P2 ](x';x") i 

V k Pl(TlP2P3 (x"; d'ld'ld'i 



(a; ; x ) 



(x ;x ) 



(3.96) 



And the diagram of Fig. 4c has ghosts on the outer legs and takes the form: 
,.#-1 - -fl-i , 10 



-« / du' / dV / du" / d 3 r" £ 



i,j,k=l 



V^^(x;d 1 ,d 2 ,d 3 ) i 
V^ a *(x'; d[,d' 2 ,d' 3 ) 



(x ; x') z ^ 3 A / r, 1 



(x ; x") 



(x';x") 



a 3 Ap 3 



(x';x") if^'V; 



(3.9c) 



The meaning of the various derivatives is determined by the subscripts - - indicating 
which propagator is being acted upon — and by the number of primes — telling with 
respect to which coordinate the derivative is being taken. For example, the derivative & 2 
of Vp^fo^ix'^d^d'^d'z) in (3.9a) differentiates the propagator i[ a2 p 2 A p2(T2 }(x';x") 
with respect to x . The overall factor of k in each expression is the one appearing in 
our definition (2.19) of the amputated 1-point function. The sign of (3.9a) is positive 
because each vertex comes with a factor of +i and there is an additional factor of +i from 
amputation: 



x 



(Mi » ) ■ ( 3 - 10 ) 

Expressions (3.9b) and (3.9c) are negative on account of the same factors with the Fermi 
minus sign of their single ghost loops. 

The various 3-3-3 diagrams were computed in six steps using the symbolic manip- 
ulation program Mathematica [12]. For diagrams (4a) and (4b) we first performed the 
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tensor algebra and acted the derivatives of the inner loop, then did the tensor algebra and 
derivatives of the outer loop and, finally, performed the integrations. The sums in diagram 
(4c) were short enough that we could do the contractions and act the derivatives for both 
loops at the same time. The tensor algebra was performed using FeynCalc [13] loaded in 
Mathematical 

3.1 Inner Loop Tensor Algebra and Derivatives. 

The inner loop is defined by the propagators which link x ,fJl and X , and by the inner 
loop derivatives $2_3 an< ^ ^2-3 • There are eight inner loop indices to contract in diagram 
(4a) and four in (4b), and there can be anywhere from zero to four inner loop derivatives 
to act. The outer derivatives, d[ and d" are retained as free vector operators at this stage. 

What one gets after contracting the inner loop indices and acting the inner loop deriva- 
tives is a series of scalar functions times the 4-index objects which we denote by the symbol, 
j-j-aipiPifi ^ ^hese are the 79 independent 4-index objects, with the appropriate symme- 
tries, which can be formed from r\^ v ', y^, t^, and up to one of each outer derivative. It is 
useful to categorize the scalar functions by which, if any, contracted outer derivatives they 
possess. What remains are functions of it', u" and y 2 . We call these coefficient functions, 
C m A(u', it", y 2 ), where the index m stands for which of the ten contracted external deriva- 
tives multiplies the function and the index A stands for which of the 79 4-index objects 
the combination multiplies. We therefore reach expressions of the following form for both 
(4a) and (4b): 

rH -i . rH -i f 43 

k / du' dV du" I d 3 r" J2 Vf^ 2U2 ^ U3 (x;d 1 ,d 2 ,d 3 ) 

Ju J Ju J 



79 



(x;x r ) i f j i3 v 3 Ap 1( j 1 (x;x") ^ Ii^ lPl 1 

A=l 



* It is not very efficient to use FeynCalc this way because most of the 2 megabytes memory of the program 
goes to define operations we never need. After completing the calculation we wrote a very short contraction 
program in collaboration with G. G. Huey which reproduces the tiny portion of FeynCalc we actually used. 
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C 1A + t- d[ C 2A + y d[ C 3A + t- d'{ C 4A + y d'{ C hA 
+ (t ■ d[) (t ■ d'() C 6A + (t ■ d[) (y ■ d'l) C 7A + (y ■ d[) (t ■ d'{) C 8A 
+ (yd[) (y d'{) C 9A + d[ ■ d'{ C 1QA 



(3.11) 



What we call "the inner loop" part begins with the sum over A. 



As an example, consider the contribution to the inner loop of (4a) from the vertex 
j = k = 41 of Table 3: 



(inner) 4 i,4i = { V^ a ^ a ^(x'; d[,d' 2 ,d' 3 ) i [ a2 p 2 A p2(j2 ] (*'; x") 



{x';x") 



= k 2 n' 2 //2 d[ ■ d' 3 d'{ ■ d'{ v ai ^ a2 v ^ a3 rfr^ T] pl ^ p2 rf 2 )^ 3 rf^ ax 
{ iA 3N 2 Va 3 (p 3 Va^ps ~ ^3*3 



-iA 



3L 



2? la 3 (p 3 Va 3 )(3 3 - 2 Va 3 /3 3 Vp 3 a 3 ] } • 



(3.12) 



After the contractions we obtain: 



(inner) 4M1 = k 2 fi' 2 fi" 2 d[ ■ d' 3 d'{ ■ d% !^A 2N iA 3N [277^ rf^ + 2 V ai ^ rf 1 ^ 1 ] 
+ (i^ 2 N «A 3L + iA 2L iA 3N ^ 

+ iA 2L iA 3L h v ai131 V piai - 3?f ^ rf^ + 5 V ai N p H ai 



_3 ?? ai/3i ^Pivi +7 j a i(pi jfijPi 



At this stage we recognize six of the 4-index objects: 

r ai/3ipicri _ ^aiPi ^pi<Ti Yl® 1 P 1 P 1(Tl = f ai t^ 1 r]P iai • 

^a 1 (3 1 p 1 a 1 _ t ( ai ^/3i)(pi ^cti) . 



(3.13) 



YjaiPiPiVi = ^Q!i(pi ^cti)/3i 



(3.14a) 
(3.146) 
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jjai/?ipi<ri _ ^ai ^pi ^<ti 
13 



(3.14c) 



There does not seem to be any point in giving the entire list of the 79 objects since we will 
not use them further. Dropping the indices and acting the derivatives gives: 

.2 



(inner) 4 i ;4 i = 



K 



26yr 4 



32x-d[ x-d'j %& x -d'{ 8x-d[ t-dj _ 8t-d[ x-dj 4t-d[ t-d'{ 



u'x ti 



2Ui + 2n 2 



+ 



( \<6x-d' x x-d'j Ad' r d'{ Ax-d[ t-d'j _ At-d[ x-d'j 2t-d[ tj% \ , , u 2 2\ 
\ u'u"x« + WW + u>u" 2 x± u' 2 u"x± + u' 2 u" 2 x 2 J [ X > 



+ 



Adj-d'j 8x-d[ x-d'j 
' u'u"x A ~ l u'u"x ti 



-3Ui + n 2 - 3n 3 - 3n 4 + 4n 9 



/ 2d[-d'j . Ax-d[ x-d'{ \ , ,„2 2n 

+ ,2 „-2 2 + ,2 ,,2 4 M H X ) 

\ w u" x l w u" a; 4 / 



5ni - 3n 2 + 5n 3 + 5n 4 - 6n 9 + 2n i3 



from which we recognize 27 non-zero coefficient functions. A few examples are: 
C 61 («', «", y 2 ) = - 6 4-^ + m jTf ] - 201n(g2 ^ 

Ul \ " ~ l7 / n lO 111111® It' It'' It® 



(3.15) 



u'W'y" 



u'u"y u 



+ 



u ,2 u" 2 y 4 



n ( l " 2\ 16 , 41n(iJ 2 y 



2„,2^ 

4 ' ' 



Cl013(« ,w ,y ) = 



#.2 Outer Loop Tensor Algebra and Derivatives. 



2„,2\ 



/ // 2n _ Aln(H 2 y 



u' 2 u" 2 y 2 



(3.16a) 

(3.166) 
(3.16c) 



What we mean by the outer loop tensor algebra is the contractions over cq, j3\; p\, 
A*2-3 an d ^2-3- What we mean by the outer loop derivatives are d[, d" and 82-3- 

The single possible derivative on the external leg acts backwards on the entire expression, 

including the limits of integration. Only its 0-component survives: 



0f — > ^ 



(3.17) 



The great advantage of the representation (3.11) we use for the inner loop is that it isolates 
the relatively simple dependence upon derivatives and indices from the very complicated 
coefficient functions C m ^(u' ,u" , y 2 ). We can contract the indices and act the various 
derivatives for each coefficient function without having to worry about its functional form. 
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An additional simplification is that we can contract the external indices, \i and u, into the 
two independent tensors — rj^ and t v — which survive the integrations. So the outer 
loop result can be expressed in terms of 2 x (10 x 79) = 1580 outer coefficient functions: 

a mA (u, u\ u", w 2 , y 2 , z 2 ; d u ) 



1mA(u, u\ u", w 2 , y 2 , z 2 ; d u ) 



The ten contracted derivatives V m {y, d[, d") are defined in expression (3.11). Because no 
derivative can appear more than once, the inner coefficient function C mA is automatically 
zero for those combinations of m and A where either d[ or d" appears in both the contracted 
derivative and the 4-index object, n^ 1 ^ 1/,lfTl (y, d[, d"). This means that only 696 of the 
coefficients (3.18) were actually used, so of course only these were computed. 

We should comment on how the outer coefficient functions (3.18) can be reduced to 
dependence upon only the three conformal times — u, u', and u" — and the three Lorentz 
squares of (3.2) — w 2 , y 2 and z 2 . First, one can re-express contractions of spatial vectors 
in terms of Lorentz contractions and 0-components: 

a-b = a-b + a°b° . (3.19) 

Second, contractions involving = tj^q can be expressed as follows: 

t-t = -l ; (3.20a) 

t ■ w = v! — u , t ■ y = u" — v! , t ■ z = u — u" . (3.206) 
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^ E vr^ U3 (x; lh,d2, d 3 ) 
i=l 

V m ^d'{) n a /^(y,d> h d'{) 

M2^2^o;i/3i] ( x: > x ') i M3^3^PiCTij ( x 'i x ") •> 



(3.18a) 
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W» E Vr^^ix; ^82,83) 



i=l 



(3.186) 



Finally, dot products between the three coordinate differences can always be expressed in 
terms of the Lorentz squares: 

1/ 2 2 , 2\ 
w-y = 2 (-iu - y + z ) , 



1/ 2 2 2\ 
y-z = ^{w -y - z ) , 

1/ 2,2 2\ 



(3.21) 



Because the outer coefficient functions involve a sum over the 43 outer vertex operators, 
they often have lengthy expressions. One of reasonable size is 7113: 

„ 119 f,. ,/ .// „,2 ..2 _2.o \ _«tf 2 * f lguV; i 8(u'-u)u'u" 8u' (u"-u)u" \ 



+ 26 



y2 r 



3mV 



10(m'-u)uV' 22m , (m ,/ -m)m / ' 



UW 4 Z 2 



24(m'-m)mV / -mK 6 M V'(-w 2 +y 2 -z 2 ) 



} ■ (3-22) 



Terms containing single factors of ln(i7 2 w 2 ) and/or ln(i7 2 z 2 ) can also appear but in these 
cases the outer loop coefficient functions are too lengthy to display. Note that the external 
leg derivative d u is really still a free operator at this stage since it can also act on the 
w-dependence in the limits of the conformal time integrations. 

The i = 34 contribution to 0:51 is in some ways more representative of what is involved 
in computing an outer loop coefficient function: 

a 5 i=V^V 3 ^^(x; 8,^2,8,) V 5 (y,d[,d'{) nf^(^,df) 



M2^2^ai/3i 



(x; x) i 



M3 ^3^10-1 



(x; x ) 



= - 2k O 2 t/ 2 ^ 3 rf*>* 8 U t-d 3 y d'{ rf 1 ^ tf^ 

- i^2L [2 r/ M2 ( ai fj Pl>2 - 2 fj aiPl fj aiPl j } 
|iA 3iV 2 ?7 M3(pi n a{)V3 - VpitTi 
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iA 



- 2k O 2 d u £ • <9 3 y- d" |l6 iA 2N iA 3N 

- 48 (iA 2N iA 3L + iA 2L zA 3iV ) + 48 iA 2L zA 3L } 



(3.23a) 



(3.236) 



kH 2 



256u'(u"~u') 256u'(u"-u) 128mV' (w 2 -?/ 2 -z 2 ) 



+ 



384 M 'K^)(w 2 -y 2 -z 2 ) 512M , K-M)u // (w 2 -y 2 -z 2 ) 



+ 6\ 



384(u"-u') 384u'(u"-u) 192u"(w 2 -t/ 2 -z 2 ) 



+ 



+ 



I92(u"-u)(w 2 -y 2 -z 2 ) 768(u" —u)u" (w 2 —y 2 — z 2 ) 



+ 



uz v 



\n(H 2 w 2 ) 



(3.23c) 



Note again that d u is a free operator at this stage. It will eventually act on all it's, including 
those in the limits of integration. 

3.3 Re- Organization of the Integrals. 

By adding to a generic integrand I(x,x',x") its reflection under x ,fl <-> x' , we can 
make the integrand symmetric. We can then enforce a canonical time ordering, u < u' < 
u" < H^ 1 , by splitting the conformal time integrations into halves with v! before and after 
u" , and changing variables: 

rH- 1 r r H- x (■ 

/ du' / dV / du" / d\" I{x, x', x") 

Ju J Ju J 

rH- 1 r pH- 1 r 

\J du' J dV J du" Jd\" [I(x,x',x") + I(x,x",x')\ (3.24a) 

I(x,x',x") +I(x,x",x')} . (3.246) 



u J Ju 

II : r r II 1 t 

du d^r' / du" I d^r" 

u J Ju 1 



The advantage of ordering the conformal times is that the three d's can be written without 
using absolute value symbols. With the ordering chosen, they are as given on Table 1. 

The generic integrands for each 3-3-3 diagram can be obtained by contracting the 
integrands of expressions (3.9a-c) alternately with rj^u (to give the a-contraction) and 
with tfj, t v (to give the 7-contraction) . We can express the integrands for diagrams (4a) 
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and (4b) by summing over the inner and outer loop coefficient functions: 



/ / // 2 2 2 q \ 
(it, it ,u ,w ,y ,z ;d u ) 

10 79 



*EE 

m=lA=l 



(u,u,u",w 2 ,y 2 ,z 2 ;d u ) C mA (u,u",y z ) . (3.25) 



/ // 2\ 



mA 



The analogous integrands for diagram (4c) are simple enough to compute without dis- 
tinguishing between inner and outer loop tensor algebra and derivatives. Performing the 
integrals gives each diagram's contribution to the a and 7 contractions of the amputated 
1-point function: 



(u) = I du 1 (fix 1 

u J Ju' 



du" / d 3 x" 



{ 


"a 




.7. 



I 1 11 2 2 2 o \ , 

{u,u ,u ,w ,y ,z ,d u ) + 



(u, u", u', z 2 ,y 2 ,w 2 ;d u ) } ■ (3.26) 



(Note that the free operator d u acts outside the integrals.) And the diagram's contributions 
to the coefficients a{u) and c(u) of (2.19) are: 



a(u) = I a(u) + 7(tt) 
c{u) = 7(11) . 



(3.27a) 
(3.276) 



It is useful to extract the universal constant factors which result from the initial k, the 
three vertices, and the four propagators: 



/ k \3 / H 2 \4 K ^H 2 

K X ( — I X 



\H 2 J V2W 2 12 tt 8 ' 



(3.28) 



It is also useful is to split the various integrands up according to which, if any, of the three 
possible logarithms they contain. Recall from sub-section 2.9 that at most one logarithm 
from the propagators can survive differentiation. We therefore have four possibilities: a 
single factor of ln(H 2 w 2 ), a single factor of ln(H 2 y 2 ), a single factor of ln(H 2 z 2 ), or no 
logarithms at all. A final re-organization is to break the integrands up into monomials of 
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the three conformal times and the three Lorentz squares. What results is a dauntingly 



diverse series of integrals having the general form: 



# 



k a H 2 



212 



7V C 



rH- 1 
(du) ' 1 / du' 



H- 1 



du' 



( u y (u')j (u") k 



d 3 / J d\" 



In ' 1 



H z (w 2 ,y 2 ,z 2 ^ 



(3.29) 



(w 2 ) 1 (y 2 ) m (z 2 ) n 

The exponents of the various conformal times may be negative, and are always < 2. The 
exponents of the Lorentz squares may also be negative. Those for w 2 and z 2 can range 
from +3 to —2, while the one for y 2 can range from +6 to —2. It should also be noted 
that dimensional analysis constrains the sum [i + j + k + 2(£ + m + n)] to be either eleven, 
if there is an external d u , or twelve if there is not. 

3.4 The Angular Integrations. 

When the vectors f ' and f " are written using a polar coordinate system in which the 
z axis of r " is along r the three Lorentz squares have the following simple expressions: 

:1 - r' 2 - (d w - it) 2 = (r' - d w + it) (/ + d w - ie) , (3.30a) 

(3.306) 

y ' = r '* - 2r'r" cos(8") + r"* - (d y - ie) z . (3.30c) 
Since the various integrands depend only on the f" zenith angle, we can perform the other 



w 



z 2 = r " 2 - (d z -ie) 2 = (r" - d z + ie) (r" + d z - ie) , 
,2 ./2 njn „raii\ , //2 /, • ->2 



angular integrations trivially: 

-2tt 



-2tt 



/ d(j)' / d9' sin(6' / ) / d<f)" = 8n 2 

Jo Jo Jo 



(3.31) 



This factor is universal and we multiply it into the other universal factors (3.28). When 

there is no ln(H 2 y 2 ) and m/1 the integral over 9" gives: 
3// sin(0 -1 



Jo 



dff 



(y 2 ) 



2\m 



2(m — l)r'r" 
-1 

2(m - l)r'r" 



/2 



2r'r"cos(6") + (r") 2 -(dy-ie) 2 ^ 1 ™ (3.32a) 



{r' + r") 2 -{d y -ief 



1—m 



(r' - r") 2 - (d y - ief 



1—m 



(3.326) 
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When there is a \n(H 2 y 2 ) and m^l the integral over 6" gives: 



-1 



In 



(r' + r") 2 - (d y - ie) 



+ 



2(m - l)r'r" \ \(r' + r") 2 - (d y - ie) 2 ]™- 1 
When m = 1 one gets a new logarithm: 

d9" Bin(0 1 = m[(r' + r") 2 - (d, - *e) 2 

//x 1*(#V) 



— r 



— r 



(3.33) 



— r 



—r 



^{ln 2 [(r' + r») 2 -K-«) 2 ]-(r"--r»)} 



(3.34a) 
> ' (3.346) 



/ ?r ^ // sin(0") 

JO y z Ar'r' 

In this case we always eliminate the extra logarithm by partially integrating on one of the 
radial variables or conformal times. 

It is useful at this stage to absorb the lower limit of the angular integration — the 
r" — > — r" term in (3.32b), (3.33) and (3.34a-b) — into an extension in the ranges of the 
radial integrations. For the case of a \n(H 2 y 2 ) and m^l the extension goes as follows: 



oo roc 

dr'r' 2 / dr"r" 2 



-1 



x 



{w 2 Y (z 2 ) n 2(m - l)r'r" 



In 



H 2 (r' + r") 2 - (dy-ie)' 



[(r' + r"f - (d y - ie) 2 ] 



) )\ + m^I _ / „ 
2]m-l V 



2(m 



/ dr' r' / dr" r" 



In 



H 2 [(r' + r") 2 -(dy-ie 



+ 



m—l 



(w 2 Y [(r' + r") 2 - (d y - ie) 2 ] m ~ l (z 2 ) n 



1 



dr' r 



dr" r" 



In 



H 2^ r > + r » ) 2_( dy _ ie) 2^ 



+ 



m—l 



4(m-l)y_ 0O ' ' ■ ' (w 2 f[(r' + r") 2 -(dy-ie) 2 ] m - 1 (z 2 y 

The reduction is similar for the other cases. 



(3.35a) 
(3.356) 



We denote by Y 2 the combination of terms which descends from y 2 after performing 
the angular integrations and extending the radial ranges: 



Y 2 = (r' + r") 2 - (d y - ie) 2 = (r' + r" - d y + ie) (r' + r" - dy - ie) . 
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(3.36) 



At this stage the answer is a long series of integrals of the form: 



# 



K^H 2 

29 



zt2 r H~ 



du 



H- 1 



du 



1 



( u y {u')j ( u ") k 



/oo rOC 
dr' / dr" Q 
-oo J— oo 



In 1 



0,1 



H z (w 2 ,Y 2 ,z 2 ) 



(3.37) 



(w 2 ) e (Y 2 ) m (z 2 ) n ' 

where Q is a quadratic function of r' and r" whose precise form depends upon what, if 
any, partial integrations were done to eliminate extra logarithms. 

3.5 The Radial Integrations. 

Since w 2 , z 2 and Y 2 factorize into products of linear functions of r' and r", the method 
of contours is especially effective. Note that the single possible logarithm can always be 
decomposed into a part which is analytic in the upper half plane and another part which 
is analytic in the lower half plane, for instance: 



ln(IfV) = In H{r" - d z + ie)] + In H{r" + d z - ie)] 



(3.38) 



The factor of (2tt) 2 which comes from the two contour integrations is universal and is 
multiplied in with the rest. 

As an example, consider (3.35b) with £ = n = 1 and m = 2: 

ln(U 2 y 2 ) + 1 



dr' r 



2 



7 // // 

dr r 



-oo 



w 



2 z 2 y2 



n 



4 (d w + dz + dy — Sie) (d y — ie) 



+ 



TT 2 d v ln 



H(d w + d z + dy — 3ze) — (d w + d z ) ln 2H(d y — ie 



+ 



TT 2 d V ln 



(d w + d z + d y - 3ze) (d w + d z - d y - ie) (d y - ie) 
—H(d w + d z + d y — 3ie) — (d w + d z ) ln —2H(d y — ie) 



(3.39) 



4 {dw + d z + d y - Sie) {d w + d z - d y - ie) {d y - ie) 

Note that there is no pole for d y = d w + d z ; in this case the numerator vanishes as well as 
the denominator and the residue is finite. Note also that this integral can be differentiated 
with respect to d w and d z to give a generating function for integrals with higher £ and n 
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values. In fact only a handful of the radial integrations really need to be computed directly, 
although this was done anyway as a check on accuracy. 

For large exponents the results can be quite a bit longer than (3.39). In the general 
case one also has to expect logarithms and inverse powers of d w — ie and d z — ie. The form 
we reach after doing the radial integrations is a long sum of conformal time integrations 
of the following type: 

K 4 H 2 r H- x rH- 1 i 



A n Ju Ju' 



{uY{ u y{u") k 



In' 



0,1 



±H (2d w — 2ie , 2d z — 2ie , 2d y — 2ie , d w + d z + d y — 3ie) 



(d w - ie)^ {d z - ie) x (d y - ie)^ (d w + d z + d y - 3ie) a (d w + d z - d y - ie) 1 
3.6 The Conformal Time Integrations. 



(3.40) 



It is at this stage that the "±" variations of the vertices become important. One has 
to sum (3.40) over the four possibilities, with the cfs assigned as in Table 1, and with a 
factor of —1 for each "— " vertex. One can see from Table 1 that the three cfs are always 
plus or minus the associated coordinate differences, for example, d z = ±(u' — u). However, 
the two sums which can occur behave quite differently for the "++" and " " variations: 

d w + d z + d y = ±2 (u" — u) , d w + d z — d y = ±2 (u f — u) ; (3.41a) 

than for the "H — " and " — h" variations: 

d w + d z + d y = ±2 {u" — u) , d w + d z — d y = . (3.416) 

As noted above, apparent poles in at d y = d w +d z are always spurious. They are evaluated, 
for the "H — " and " — h" variations, by first taking the limit dy — > d w + d z , and then 
substituting for d w and d z from Table 1. 

The obvious strategy at this point is to decompose the integrands by partial fractions. 
Because only single powers of logarithms arise, this will always suffice to reduce the result 
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to a single integral.* The remaining integral can also be decomposed by partial fractions, 
however, an integral of the form: 

ln(x — a) 



I 



dx- 



x — b 



(3.42) 



cannot be expressed in terms of elementary functions for a 7^ b. When this occurs the 
strategy is to first extract any ultraviolet divergences by partial integration, then expand 
the logarithm and integrate termwise. Since we only require the leading order form as 
u — > + , this is straightforward. Consider, for example, the following: 



Su 



H 1 ]r\( Hv r ) H^ 1 
du'^ '— = In(Hu') \n[H(u' -u- ie)] 



u — u — ie 



"tf" 1 du' 



= -\n(Hu) In(-iHe)- [ H ^r\HHu')-Y^Uf) n ] (3.436) 

Ju U l \ / J 

n=l 

00 

= - \n(Hu) In(-iHe) + \ \n 2 (Hu) + ( Hu "> U • ( 3 - 43c ) 

n=l 

Of course one can use these methods to perform the integrations in the original order and 
with the original variables of expression (3.40). We did it both ways and compared each 
term as a check on accuracy. 

We have selected, as an example, the term which emerges from step 3 above in the 
form: 



/ — \n[H(u' - u)} (3.43a) 

Ju 



2 // 



d\> I d\" y 



2 4 2 



It has already been shown that the integrations over r ' and r " give: 

1 ^ ^ f H 1 ^ r*- 1 du" f 1 

16 2 7 7r 4 u J u J u i u 2 u" \ (d w + dz + dy — 3ie) (d y — ie) 



(3.44) 



+ 
+ 



dy In 



H(d w + d z + d y — 3ze) — (d w + d z ) In 2H(d y — ie) 



(d w + d z + d y - 3ie) (d w + d z - d y - ie) (d y - ie) 



(3.45) 



* Although sometimes not in terms of either v! or u" . 
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If we neglect terms of order e, the "++" and "H — " variations give the following expressions: 
u 1 k A H 2 a ft 1 ' 1 J , f H ~ X du" f 1 1 



16 2 7 tt 4 Ju J u > u 2 u" 1 2 (u" - u' - ie) (u" - u - ie) 

1 ln[2if(u" - u - ie)] 1 \n[2H(u" - v! - ie)] 



4 [u" — u — ie) (u r — u — ie) 4 (u" — u — ie) (u" — u' — ie) 
1 ln[2H(u" - v! - ie)] 
4 (u" — u' — ie) (u' — u — ie) 



(H -> -if) I , (3.46a) 



-#1^4, r 1 *- r^ M » in '- 2g (";-";: fe )i, + ^r" g) ■ ^ 

16 2 7 tt 4 7 U 7 U / 2« 2 w"(M"-«' + ze) 2 v y 

The " " and " — h" variations can be obtained by complex conjugation of (3.46a) and 

(3.46b), respectively. 

The first integral in the "++" variation (3.46a) has a simple result: 

x /-if- 1 7 /-if" 1 du // 

du / 777777777 77 77777777 — 7 777 = (3.47a) 



2w 7 U J u ' u"{u" — u' — ie) (u" — u — ie) 



£J £ + i{-'" 2 (ir) + 2M-«e) ln(l - Hu) - In»(l - *«)} 



n=l 

Note how it illustrates the arguments given in sub-section 2.7 for the inevitability of infrared 
logarithms when an £-loop graph contains t logarithmic ultraviolet divergences. Nice as it 
is to have exact expressions, what we really want is an asymptotic expansion for late times 
(u^0+): 



1 f H , / f H du" In 2 (Hu) _/ln(£Tu)\ , n _ 
1 du' - 177 - Tl -. - j-j, - - = l—L + o(-^-L (3.476) 



2w 2 7 W Ju' u"(u" — u' — ie) (u" — u — ie) 4w 3 V 

The asymptotic expansion for all of (3.46a) is: 

„ k 4 H 2 f 1 ln 3 (M 1 r , . Z7T1 ln 2 (M ^ /ln(£Tu)\ 1 , , , 

The analogous expansion for the "H — " variation (3.46b) is: 
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Summing these, along with the respective complex conjugates for the "- 
variations, gives: 



# 



k a H 2 



2 7 7T 4 1 16 



1 ln 3 (Hu) „ 7r 7 3 , , n . 

+ S ^32 + 16 ln(2) 



In 2 (ffu) | /ln(ffu) 



-" and "-+" 



(3.49) 



Note the triple log terms. From the discussion of sub-section 2.8 we know that these cannot 
appear in the full result,* however, they do appear in pieces of it such as this example. An 
important check on accuracy is that the sum of all such triple log terms cancels. 

Many, many terms such as (3.44) emerge from step 3. For example, the 7 contraction 
of diagram (4a) — which is also diagram (2a) — consists of 842 terms whose integrands 
contain a factor of ln(if 2 ?/ 2 ), 301 terms which contain ln(H 2 w 2 ), another 301 terms which 
contain ln(if 2 2; 2 ), and 1663 terms which contain no logarithm at all. When all the 3-3-3 
integrals are computed and the results summed, the totals are: 

K±H 2 



a(u) = 
7(u) = 



2^4 
2V 



. ln 2 (M ^An(Hu) 
-492 + 234 + 216 V^ + Of 



■ir 



u 



(3.50a) 



1157 



400 



56 ]^ + o(!^)}. 



6 3 i w \ w 

where the three numbers between the square brackets refer to the contributions from 
diagrams (2a), (2b) and (2c), respectively. Taking account of (3.27) this gives the following 
3-3-3 contributions to the coefficient functions a(u) and c(u): 



K±H 2 

a 333(«) = 
c 333(«) = ^4 



1795 302 I6O1 
+ -r- + 



+ 



18 
1157 



9 
400 



6 



3 

- 56 



vr 



\n 2 (Hu) | ^n(Hu) 



vr 



if 



(3.51a) 
(3.516) 



The results quoted above contain the contributions of some degenerate terms whose 
reduction is actually much more akin to that of diagram (2d), to be described in the next 



* In fact the argument of 2.8 applies for the result from any triplet of vertex operators. We checked this 
explicitly for the triplet in which the x^ vertex operator is #10 and the vertex operators at x'^ and x"^ 
are both #41. Step 3 in the reduction of this triplet gives 117 terms containing a factor of \n{H 2 y 2 ) and 77 
terms with no logarithms. 
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section. These terms come from the action of two 0-component derivatives on the normal 
part (2.14a) of an outer leg propagator which connects two "+" vertices. The vertices can 
be either x^ and x ,fJl or and x' ,fX '. What such a double derivative really gives is: 



d d'o iA N (x; x') = 2 d d' 

H 2 



u u 



+ 



u u 



Ax — \ Au\ + ie Ax + \ Au\ — ie 

1 1 

+ 



(3.52a) 



8tt 2 Ax \ Ax - \Au\ + ie Ax + \Au\ - ie 

Au 



+ 



9{Au) - e(-Au) 
2uu' 



Au 



(Ax - | Ait | +ie) 2 (Ax + \Au\ - ie) 2 
2uu' 



(Ax - \Au\ + ie) 3 (Ax + \Au\ - ie) 3 



2 8(Au) 



u 



u 



(3.526) 



(Au + ie) 2 (Ax-ie) 2 _ 

(Recall that Au = u' — u and Ax = \\x ' — x\\.) The reduction procedure outlined in 
this section accounts for all but the terms proportional to S(Au). They are of order e but 
sufficiently singular at x = x ' that the result is a spacetime delta function:* 



H 2 u 2 8(Au) 



An 2 Ax 



1 



1 



(Au + ie) 2 (Ax-ie)' 



H 2 u 2 ie 5(Au) 



7T 2 (Ax 2 + e 2 ) 2 

n 2 id 4 (x'-x) . 



(3.53a) 
(3.536) 



Performing the now trivial integration over x'^ causes the 3-3-3 diagram to degenerate to 
the 4-3 topology of diagram (2d). The reduction thereafter is the same as will be described 
in the next section. Of course delta function terms can happen whenever a propagator is 
doubly differentiated. They give 4-3 diagrams when either of the two outer leg propagators 
is affected; when one of the two inner propagators is affected the degenerate diagram has 
the topology of figure (2e) and fails to contribute to leading order for the same reason. 



Similar terms of order e arise from the logarithm part of the propagator but they are not singular enough 
to survive in the unregulated limit. 
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4. The 4-3 Diagram 



The 4-3 diagram consists of an outer 4-point vertex joined to a freely integrated 3-point 
vertex as shown in Fig. 5 below. From the discussion of sub-section 2.7 we see that this 
diagram can contribute at leading order when a factor of ln(Hu) from integrating the free 
interaction vertex combines with a factor of ha(Hu) from an undifferentiated propagator 
logarithm. The physical origin of the first logarithm is the growth in the invariant volume 
of the past lightcone, while the second logarithm comes from the increasing correlation of 
the free graviton vacuum of an inflating universe. 



Fig. 5: The tensor structure of the 4-3 diagram. 

The 4-3 diagram has an obvious 3-fold symmetry, corresponding to permutations of 
lines 1, 2 and 3. In computing such a highly symmetric diagram it is not efficient to fully 
symmetrize the vertex operators and then divide by the symmetry factor of 6. The better 
strategy is to symmetrize the 4-point vertex operator only on the external line, and then 
contract this partially symmetrized vertex operator, through propagators, into the fully 
symmetrized 3-point vertex operator. 

We can read off the graviton 4-point interaction from expression (2.4): 



\i v 




Pl°l 



P3°3 
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4 4 l = ^ fi2 { W v> V" M - ^ 2 - W <w + W w 

- |v> ^ a/? V"" + ^ a/3 V" - ^ ^ a/? ^ V 

~ /a V> W + 1^3 V^ 7 ta ~ l^B ^ V,« V^ 7 



+ 



J /a V> W + ^ V^ 7 V> V>, 7 " 3^ V^ 7 <W ^ 

1^ V^ 7 V>a 7 ,M V" M " ^ V^ 7 V>a 7 ,M ^ + ^ V>a/f V> 7 / 



- \r Q ^ + v^ 7 i>aw - W* p ^ ta v 

- \r p v^ 7 v vw, 7 + lr p r u vw?, 7 V - \^ ^ v ^ V^ 7 
+ \r p r v tap* v> - \r p r v ^ to - \^ ^ v ^ ^ 

- \r p v v vv, 7 v< v - \r p r v vv,7 v + r p r v vv, 7 i>\ v 

+ \^ ^ V>a 7 ,M - ^ ^ V t v + ^V> ^ 1>a0 V h 

- vw v r v t v - y« w ^ r v u + ^ r v u 

+ ^ 1> ^ - \^ ^aP ^ ^ U} • (4.1) 

(Recall that ip = ip^ and = tj^q.) Almost all of this ungainly expression can be checked 
against published results [11] by taking the flat space limit: u = H~ l — t and H — > 0. The 
procedure for partially symmetrizing a vertex was described at the beginning of Section 3. 
As an example consider the final term in (4.1): 

-A 2 i; ap ^ V V K ■ (4.2) 

If the indices fx and v and the derivative d u represent the distinguished line, a valid partial 
symmetrization of this interaction gives the following vertex operators: 
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- K 2 n 2 1 d!f v u){a3 v ai{a2 v^ 1 , (4.36) 
- n 2 n 2 1 d u v ^ n^t^ V 32 )^ 3 , (4.3c) 

-K 2 n 2 I & t? 1 ^ 1 8^ r, a ^ rfr^ . (4.3d) 

These are vertex operators i = 127, 130 in Table 4. The fully symmetrized cubic vertex 
operators are given in Table 5. They were obtained by interchanging legs 2 and 3 from the 
Table 3, but we have taken account of symmetries to reduce the 2 x 43 = 86 terms to only 
75. 

The unprimed vertex at x^ is "+" type and we must sum the primed vertex over both 
"+" and "— " variations. Since the (++) and (H — ) propagators are identical for u' < u, the 
relative minus sign between the two vertex assignments allows us to restrict the range of 
integration to u < u' < H~^. In this region the (H — ) propagator is the complex conjugate 
of the (++) one, so we can write the total contribution of the 4-3 diagram as twice the 
real part of the (++) term: 

Tjg = a 43 (u) rT + CA3(u) & t v 

( rH- 1 f 130 
= 2Re|-m J du' J d 3 x' V^^^ (x; d u , d x , 6> 2 , <9 3 ) 

i (x;x') i a2 p 2 Ap 20 - 2 J (x; x') i a3 ^ 3 Ap 3<73 (x;x') 

75 . 

£ Vf^^V^A^) . (4.4) 
3=1 J 

The subscripts i and j refer to Tables 4 and 5 respectively. Recall that where each derivative 
acts is indicated by primes and subscripts. For example, the derivative d% in the 4-point 
vertex operator acts on the first argument of the propagator 2 [3 A3] (2;; x'). The derivative 
d u acts on all -u's in the vertex and the three propagators. 

The entire calculation was performed by computer using Mathematica [12] and Feyn- 
Calc [13]. The first step was to contract each pair of vertex operators into the three internal 
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propagator and write the results onto a file. The next step was to act the internal deriva- 
tives (<9i_3 and d[_%) and store the results for each pair of vertex operators. Selected 
vertex pairs were computed by hand to check the procedure. 



% 


Vertex Operator 




Vertex Operator 


1 


aipi g( a 2 r7 /3 2 )(a 3 M 
8 ' ' 3 ' 2 


66 


4t<9 7/ ^ v ri^l 31 n 012 ^ 2 r?" 3 ^ 3 1 ■ a 


2 


IB B^ n zy )( a 3 ffe) rPt-\&\ rfx-2^2 
^3 / ' / 


67 


_l r? M(«l 77/^1)^ ^2/32 -^3/53 ^ . 

§'/ '/ / / ^0 


3 


1 M (oi Mu q(<*2 /3 2 )(a 3 Q fo) 


68 


-la, r?^^ r? ai ( a2 r?^ 2 )^ 1 r?" 3 ^ 3 t ■ B-x 


4 


-ia, aS M tP$ W" 2 


69 


„M)(«i |,/3i)(a2 ?? /3 2 )(^ ^asft ^ 9 . ^„ 


5 


„//)(«! aft) a("3 /3 3 )(a 2 -jft)^ 


70 


ia, r^X"! r,0i)(<X2 Jh){v n a 3 f3 3 + . a 
2"/ / / / 


6 


a(^ n v)(ai „/9i)(a3 riftOfo^ 
^2 ' ' ' 3 


71 


la, r?^^ n a i)( a 2 r//3 2 )(a3 n (3 3 )(Pi + . q„ 


7 


9„ r>^)( a 2 77 /32)(Q!i #ft) ^(a 3 

^ Li '/ '/ ^-4 v '/ 


72 


iB^ n v ^ ai B^ 77 a 2/3 2 v a 3 p 3 
2 2 ' 3 ' ' 


8 


ft, 9^ ?7^)( a 3 M(<*1 r>ft)(o 2 £/3 2 ) 


73 


1 q iiv t (ai W?i)(a 2 ^2) a3/ 3 3 


9 


„A»)(ai r>ft)(o 3 o(a 2 (3 2 ){y 
i '/ w 2 3 ' 


74 


— \n^ v n ai ( a2 7102)01 ^303 f)~ . f)o 
4 / / / / ^ 


10 


f) 8^ n u )( a2 r^K^i rA)( a 3 fPs) 
3 ' * * 


75 


_ 1 ^(tt 2 ^2)^ yjQil/Jl n ^303 . Bo 

4 / / / / ^ 


11 


2 / ^3 ' ' 2 


76 


-ia, r7^( a2 r?' 32 ^ ?7 ai/31 ?? a3 ^ 3 t ■ a 


12 


-id^ r?^" 3 d^ 2) v aiPl 
2 2 ' ' 3 ' 


77 


-k, 77^^ Tl"!/ 3 ! ?? Q!2 ( Q! 3 Mfo t ■ B-X 

4 t* / / / / O 


13 


-lid,, 7?^)( a 2 t(a 3 ^aift 


78 


8 ' 2 3 ' ' 


14 


— IB B^ r) v )( a 3 n/ 3 3)(a 2 +ft>) ^aift 


79 


— XB^ B 1 ^ n a i/5i 7,0:2/32 7,0:3/33 
8 2 3 ' ' ' 


15 


4 / 3 ' ' 2 


80 


A u> l I I < 


16 


_1^)(« 2 ^2) ^ 9 («3 ^ 


81 


l^(a 2 ^ a (oi ^1) r/ a 3 /3 3 


17 


-\B U B^ T]^ a2 r 7 / 3 2)(a 3 t /3 3 ) v <xi(3i 


82 


1 a v ) 71 a 1 (a 2 „/3 2 )/3i ^03^3 
5 a 2 a 3 " " '' 


18 


qO 1 ^){a 3 ^/3 3 )(ai ^/3i)(a 2 ^ft) 


83 


lQ uT] fj,(ai ^{)v ^(02 ^/3 2 ) ^a 3 /3 3 


19 


^(02 ^ftXai ^1) ^("3 r7 /3 3 )( i y 


84 


^ u 77^^ a^" 1 t^ 1 ) r;02(o 3 ^s)/*! 


20 


a u ^i)(a 3 t ft) 9 («2 ?? /3 2 )( I , 


85 


1^(03 j^z)v ^axia.2 T/ /3 2 )/3i ^ . ^ 


21 


a u 0^ ^^("2 jjfoKas v (3 3 ){ai t (3i) 


86 


1^ 7/ M(o 2 ^2)1/ ^1(0-3 v P 3 )fc t . g 3 


22 


q{v ^v){a 3 ^/3 3 )(a 2 ^/3 2 )(ai #ft) 


87 


_\^iv ^aiPi ^a 2 (a 3 ^/3 3 )/3 2 ^ . 



Table 4: The partially symmetrized quartic pseudo-graviton vertex operators v r / II/Qa/3lQ2/32Q3/33 without 
the factor of k?Q? . 
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% 


Vertex Operator 


i 


Vertex Operator 


23 


Q u ^)(a 2 gfo) f(ai ^0i)(a 3 ^0 3 ){v 


88 


-jgd u i]^ a3 rp^ v Tj ai 1 r] a 202 1 . q 3 


24 


-\r]^ y rf^ rf 2 ^ 2 d^ 3 


89 


l^(ai ^0\)v ^a 2 (a 3 ^0 3 )0 2 Q 2 . 


25 


-^d u rf v rf^ rf 2 $ 2 d { 3 a3 


90 


lQ u v ^(a 3 v 3 )is ^(a, v 2 )0i t . Q 3 


26 


-j§d u rf^ 1 rf 2 ^ 2 rj^ 3 


91 


-rj^){ai v 0i)(a 3 v 03)(&2 ^AOO q 2 . g 3 


27 


^0i)v ^a 2 02 Q^ 3 Q0z) 


92 


—d u 7?^X"2 T j02)(a 1 ^0i){a 3 ^0 3 ){v t . 


28 


\d u rf v ?7 a i( a 2 r] 02)0\ q^ 3 ^3) 


93 


— \d^ rf^ 011 r7 a2 ( a 3 ^03)02 


29 


?7 ai ( a2 77W1 77^ 


94 


-\d u J}^ 3 v 3 )v t (ai v 0i)(a 2 d 02) 


30 


^0(02 jfr){y Q^ 3 ^3) 


95 


l^iv ^a 1 (a 2 r j0 2 ){a 3 ^03){01 Q 2 . Q 3 


31 


~2^u r]^ ai r]^ 012 rP 2 ^ v d^ 3 t^) 


96 


l^)(a 2 r j02)(a 3 ^0 3 ){v ^ai0i q 2 . 


32 


-\d u t y ) n ai ^ a2 r]^ 013 77&)(/?i 


97 


\du 77^X^2 Tj 2 ){a 3 ^03){v ^101 t . g 3 


33 


^) (ai ^)(a 3 fi^) v a 2 ? 2 


98 


irT dt 4 l] V a2{a3 ^ 2 


34 




99 


\d^2 d V 3 r]® 1 ^ 1 rj a2 ( a3 rjP 3 ^ 2 


35 


-\d u tf v t( ai n^ a2 77^2 )( a 3 


100 


\d u i]^ a3 rp 3 ^ 11 dg" 2 t^ 2 ) ^101 


36 




101 


_l r] fi)(a 2 7? /3 2 )(ai ^0\){ol3 ^03)^ Q 2 . Q 3 


37 


_l^)(ai } d^ 2 ?/2)^ 77 a 3/3 3 


102 


-\d u T?^ 1 )^ 3 v 03)(a 2 v 02){v t . Q 3 


38 


^^)(ai ^i)(a 2 ^a 3 /3 3 


103 


_l r7 M)("2 v 2 )(a 3 v 3 )(u q0i) 


39 


^7 MZy ^3° 2 t/^X"! 77ft )( a 3 ^3) 


104 


-\d^2 d v 3 77 a l)( a 2 ^02){OL 3 ^0 3 ){0l 


40 


7^)02 r7 /3 2 )(a 3 ^3) ^ai/3i 


105 


-d u r]^ 012 r^Xas ^s)^ g( ai f0i) 


41 


J77^ 77 ai ( a2 V P2)P1 9 («3 ^3) 


106 




42 


1^(02 yfi)v v aipi 9 («3 ^3) 


107 


j^V-v ^a\0\ ^a.202 fa 3 ^fi 3 


43 


^(a 2 ^/3 2 )^ ^aiA ^( a 3 ^ 3 ) 


108 


— J_ 7(7* f) v } r> a ^0i n a 202 r,O! 3 3 


44 


\d u rf^ 1 r ? a2 ( Q! 3 ^03)02 


109 


^.^(ai JJ01)" yf-202 Q^ 3 t 03) 



Table 4: {continued from previous page) 



47 



i 


Vertex Operator 


i 


Vertex Operator 


45 


\r]^ v d^ 1 ^)( a 3 qPz> ^202 


110 


1 rgUV r jai(a2 ^02)01 fCt3 £03 


46 


i] l/ ^ a3 d^ r] ai ^ 1 r} a2 @ 2 


111 


■g-t^ d^ rj ai ( a2 rp 2 }^ 1 rj a3 @ 3 


47 




112 


-2.rf)(ai v Pi)(a 2 v 02)(v Q^ 3 t 0s) 


48 


\d u t& rf^ a2 df ] v ai/31 V a3p3 


113 


-±d u T]^ ai ^i)( a 2 jfr)^ f a 3 t fo 


49 


\r\V" v dc, ai rjP 1 ^® 3 r\ a2 ^ 2 


114 


T ja 1 )(a 2 r jP2)(a 3 r j03){0i 


50 


rf^ 013 d^ r\ ai ^ x rj a2 ^ 2 


115 


—-dlf 77") ( ai t/^ 1 )^ 3 t^ 3 ) ^202 


51 




116 


_l 7] fi)(ai qPi) f (a 3 v a 2 (32 


52 


\d u d^ 7]^ a2 th) rf^ rf 3 ^ 3 


117 




53 




118 


— l^G" ^)( a i ^i)(a 2 ^2) ^a 3 /? 3 


54 


-\d^ rf^ 013 d^ rf 1 ^ 012 rf 2 ^ 1 


119 


^ai(a 2 ^/? 2 )ai Q^ 3 f0 3 ) 


55 


-\d u T]^ ai t^ 2 7]^ a3 


120 


1_ V K*2 rfr)u rfifh d ^ 3 tfh) 


56 


-\d u rf^ a2 d^ rj ai ^ a3 r]^ 3 ^ 1 


121 


Tj-du 7/^( a2 ^2)^ ^ai/?i ^3 ^3 


57 




122 


1 r) V ^ n a l/3l n a 2(a3 n 03)02 
2w 3 ' ' ' 


58 


_ 1^ ^)( a 3 ^ai(«2 ^02)01 


123 


^-77/^ ^( ai ^Pl){a3 f03.) ^202 


59 


-\d u v (3i> d { * 2 rfiK a * 


124 


^-3^ r] v ^ a3 r] aif31 rf 2 ^ 2 


60 


-\d u d^ ri^( a3 rf^ a2 rf 2 ^ 


125 


^-d u ff ,v rj 011 ^ 1 t( a2 77^ 2 )( a 3 ^3) 


61 


— \d^2 rf^ 013 ^ 3 )( ai d^\ a2 ^ 2 


126 


1 +(n r,v)(a2 p)@ 2 ) rp-X&X ^0.3^3 
2u ' 3 ' ' 


62 


— \d^ 77") ( tt3 ^ 3 )( ai i] a2 ^ 2 


127 




63 


-\d u 7]^ rjt 3 ^ 3 rfrK a *d^ 


128 


^^)("3 ^3) ^«i(a 2 ^/? 2 )/3i 


64 


-\d u r}^ ai d^ 2 rj^ h 'r) a3 $ 3 


129 


^1)^ ^("2 ^2){a3 f(3 3 ) 


65 


^•yf v fqUlfil ^202 ^303 q 2 . 


130 


_It(A* fj^^l ^"2(0:3 ^03)02 



Table 4: (continued from previous page) 
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3 


Vertex Operator 


J 


Vertex Operator 


1 


-i?v Piai v P2<72 d!f 3 t a ^ 


39 


_l^P2CT2 ^(Pl ^eri)(7 3 ^ . 


2 




40 


T jPl){P2 ^2){P3 ^03) (0"l ^ . ^ 


3 


_ 1 V P3<T3 V P1<?1 t °2) 


41 


^Pl)(P2 ^XpS ^CT 3 )(o-i ^ . ^ 


4 


^■qPi{P2 ^ct 2 )cti q(j>3 t <7 3 ) 


42 




5 




43 


\d { z P2 dl 2) ^ (PI jfn)0Z 


6 


ti' 1 


44 


_ 1 „Plffl „P3^3 ^ t^) 

Zt£ O 


7 


-l^PS „CT3)(Pl ^) „P20"2 
1£ Z 


45 


_ 1 „P3^3 „P2^2 ^(Pl ^l) 

Zt£ Z 


8 


-l^Pl „^l)(P2 „P3<T 3 


46 


_ 1 „P2^2 fl(P3 t <7 3 ) 

Zt£ 1 


9 




47 


Ij^PliPi ^(P2 ^cr 2 ) 


10 


1^(71 ^2 ??( T 2 )(p3 9 J3) 


48 


l J ^P3{p2 ^2)03 giP 1 fVl) 


11 


1^2^ a (P3 7? a 3 )(pi Q*l) 


49 


1^2 (pi „cri)cr2 ^ £<7 3 ) 


12 


lrjP3<T3 d[ Pl <l)(P2 af 2) 


50 


r? " 1 ^ 2 t CJ 2) r? P30-3 


13 


-d{ Pl 71 CTl )(P2 r? CT2)(P3 

O ' ' Z 


51 


U Z ' ' 


14 


-a! P2 r7 CT 2)(P3 „ct 3 )(pi 0^1) 


52 


ljd\ P2 77 <72 )(P3 picn 
11 1 ' ' 


15 


_q(P3 v a 3 )( Pl „cti)(p2 ^f 2 ) 
2 ' ' 1 


53 


_ a (pi „<7i)(p 3 „a 3 )(p 2 ^2) 

Z ' ' 


16 


_^(P2 <T 2 )(pi „(T1)(P3 d^3) 
O ' ' z 


54 


_^(P3 a 3 )(p 2 ^Xpi 9^1) 

1 ' ' z 


17 


_^(P3 ^7 3 )(P2 ^cr 2 )(pi ^l) 


55 


_^(P2 ^ 2 )(pi ^i)(p 3 


18 


-<9^ Pl 77 CT l)(/>3 ^CT3)(P2 ^2) 


56 


-l^Pl^l ^3 flO* 3*2) 


19 


-l^i^i ^2 a (P3 ^3) 


57 


-1^3 ^2*2 



Table 5: The fully symmetrized cubic pseudo-graviton vertex operators yP iai P 2a2 P scr3 without 
the factor of nQ' 2 . 
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3 


Vertex Operator 


3 


Vertex Operator 


20 


rfc * * O 1 


58 


_ i„P2<r 2 „pi^i 4 P3 af 3) 

r± * * 1 


21 




59 


l^l(P3 ^ 3 )ai a (P2 ^2) 


22 




60 


1^3 (P2 r7 <72)<73 fljO 


23 


1 77 P2(P3 ^3^2 ^(Pl 


61 


i^2(pi ^1)^2 ^ a° 3) 


24 


l v Ps(j>i d[ P2 d°f 


62 


l^ 1 ^i)(P2 9^2) ^3 


25 




63 




26 


1 5 (P2 ^a)^! 0*0 ^3 


64 




27 


1^3 ^3)^2 ^2) v p 1(7l 


65 


1^3 ??( T3)(P1 fljO V P2CT2 


28 


Iq(Pi „<71)(P2 p 3 a 3 
2 2' J ' 


66 


l^ 2 ^2)(P3 ^3) Piai 
2 1' 2 ' 


29 
— ■ / 


lfl(P2 <7 2 )(p3 flfs) p 1( Jl 


67 


i(9l Pl r? cri )(P2 #f 2 ) riPscrs 


Qfl 

oU 




DO 


1<kiPiCi ^,02(72 r.P30~3 f) f) 


«J J. 


InPiVi n p2^2 n P3<T3 a . a 
5'/ '/ '/ u 2 u i 


69 


-l n Pl(P3 n CT3)0"i „P202 .% . A, 

5'/ '/ '/ ^2 ^3 


32 


l n Pl<Tl n P2CT2 nPZO-i f) . f) 


70 


_l r7 P3(P2 71 0-2)cr 3 ^piCTl Q . 


33 


2 ' i » z o 


71 


_l r? Pl(P2 77 <72)o-l „p 3 0- 3 « . « 


34 


-l n P2(p3 r pz)^2 n Pl<Tl Q„ . Q, 
2 f 1 1 o i 


72 


_1„P1^1 „P2^2 «(P3 ^3) 
2 ' ' 12 


35 


1^3 (Pi „Cri)<7- 3 p 2 (T2 Q . Q 


73 


_l^Pl(P2 ^2)^1 ^303 ^ . ^ 


36 




74 


^Pl)(P2 ^2)(P3 ^3)^1 Q x . Q 2 


37 


-\ d i P2 d \ 2) v p3a3 v piai 


75 




38 


_ l^Pi^i ^P2(P3 ^ct 3 )<72 g 2 . 







Table 5: (continued from previous page) 
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As an example, consider the contraction of the i = 1 and j = 10 vertex operators: 



i« ^^^(x; c\, d h d 2 , da) i L ft A pi J (x; a/) • LfcA^J (x; x') 



(x;x') y/T P2W V;<9i,<9^) 



"TO 



i A lN 2 V ai ( Pl V ai )p! ~ Varf! Vpiax ~i&lL 2 V ai ( Pl - Va^ V Pl a, 



iA 



2N 



2r la 2 (p 2 Va 2 )p 2 - ^ 2 /3 2 Vp2<T2 - iA 2L 2 Va 2 (p 2 Va 2 (i 2 ~ 2 Va 2 {3 2 Vp 2 a 2 



iA 



3N 



iA 



3L 



2l la 3 (p 3 Va 3 (3 3 - 2 Va 3 {3 3 Vp 3 a 3 \ (4-5a) 



2r M/> 3 TIadfo - Va 3 (3 3 V P3 * 3 
-f k a fi 2 fi /2 rf v (-2 iA 1N + 3 iA 1L ) 

49 2 • <9 2 <9 3 • <9 3 - 2d 2 • • d 2 + 2c>2 • d 3 d' 2 ■ d 3 ] iA 2N iA m 

+ 

-2d 2 -d 2 d 3 - d' 3 + 3d 2 ■ d' 3 d 3 ■ d' 2 - 3d 2 ■ d 3 d 2 ■ <9 3 iA 2L iA 3N 



+ 
+ 

2^u 2 u f2 



2d 2 • <9 2 <9 3 • d 3 + 3<9 2 ■d i d 3 -d 2 - 3<9 2 ■ <9 3 d 2 ■ <9 3 iA 2N iA iL 



d 2 • <9 2 d 3 • d 3 - Ad 2 • <9 3 d 3 • d 2 + 5<9 2 • <9 3 d 2 ■ d 3 iA 2L iA 3L \ (4.56) 



T] 



px> 



24^- + 24 

w 8 



u 2 t-w 2 

w 8 



- 336 



uv! t-w 2 

w 8 



+ 24 



u' 2 t-w 2 AO U t-W r> A u 2 



+ 48^-^ - 96^ - 24^C + 24^ 
w° w° w b w 4 



2„..2< 



- 48^^ - 24^ 
w b w° 



ln(H z w z ) 



+ 



uu 1 t-w ,<,uV t-w 



"32 T77 

w w 



32^ 



2 /2 . 2 /■> t 

U U t-W nnUU t-W 



-to h 448 nr 

w w w w 



32- 



3 



-jo h 64 g- 

yjW yjS 



u 2 u' t-W 



+ 32^t - 64 "" /2 8 f,MJ + 128^ + 32^C - 32^ 

w s w s w s w s w° 



(4.5c) 



Note that we have dropped the order e terms in derivatives of w A 



8 2 

w = 2 w,i — 2i e t. 



2 w 



(4.6) 



The neglected terms affect ultraviolet divergences but not the infrared terms of interest. 



The next step is to perform the spatial integrations. The angular integrations give 
a factor of An and the radial integral was done by the method of contours. When no 
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logarithm is present we obtain: 



/ | POO 



(r — wq + ie) N (r + wq — ie) 



N 



, (-1)^-1 (2iV-5)!! / 1 
71 % 2 N + X (N - 1)! \w - ie) 



(4.7a) 



The logarithm gives: 



,2JV 



8tt 2 z 



/ 



r 2 In H 2 (r — t^o + ie) (r + — ie) 
(r — wq + ie) N (r + wq — ie) N 



+ 



2 N+l (AT _!)! 
~ V 



| In ^2H(w - ie) 



+ In 



-2H(w Q - ie)} } ( -) 

J > \wq — ie/ 



(4.76) 



1 \27V-3 



8n 2 i(-l) N ^} (2N-2- k)\ [2N-2- k- k 2 ] / 1 



— 3 — k) \wn — ie) 



2 27V-1 (iv - l)! ^ k (N - 1 - fc)! (2N-2- k)(2N - 3 - k) \w - ie, 

K — 1 

The case of N = 2 has to be treated specially: 



/ (fix 1 — A = -8tt 2 z - — 

J to 4 8 tog — 



ie 



3 / hi(# 2 «; 2 ) 



w 1 



-a* - 



ln[2H(w -ie)} + ln[-2g(w -ie)1 



wo—ie 



+ 



WQ-ie J 



(4.8a) 
(4- 



Performing the spatial integrations for our i = 1, j = 10 vertex pair gives: 



• 0C J 



a 3 /3 3 A P30-3 



ft4ij " 2 nv j M2H(u'-u-ie)] + ln[-2H{u'-u-ie)] 7 1 



2^u 2 u' 2 



u'—u—ie 



2" u'—u—ie 



, 15 u' m _ 7 it' 2 M 2 v 
"•" -g- ( u /_ u _, e )3 (u'-u-ie) 6 f • 



(4.9) 



The final step is to integrate over the conformal time of the free vertex. The integrand 
always consists of a numerator — which may contain a pair of logarithms — and a de- 
nominator — which contains up to two powers of u , and up to seven powers of u — u — ie. 
Decomposing the denominator by partial fractions results in four distinct terms requiring 
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special attention: (u' — u — ie), u', u' 2 and (u' — u — ie) k for k > 1. When the pair of 



logarithms is present the four denominators give 



(i) 



/ W=£=Te {^[2H(u' -u- ie)] + ln[-2#(u' -u- ie)] } (4.10) 



= ^jln 2 ^^' 



u — it) 



-2H(u — u — ie | 



In 2 [2(1 - ifu)] - ln 2 (2#e) + in In 2(1 - Hu) - \n 2 , 



(ii) 



u — ie 



{in 2#(u' - u - ie)] + ln[-2#(?/ 

00 , i^- 1 

> {ln 2 (2Hu , ) + inHHu , )+2j2l?(%) } 

fc=l u 

-2 ln(2) ln(ifu) - ln 2 (ifw) - in \n(Hu) + 2 [Hu$(Hu, 2, 1) - C(2) 

-2H(u - u-ie) } 



(4.11) 



+ ln 



(hi) ^ ^{\n[2H(u'-u-ie) 

— > (-4 ln(2W) - % + 2^ ln[l - ^1) 
= lln(^) + g(l-^){|7T + ln[2(l-^)]} , 



(4.12) 



(iv) / 

Ju 



du 1 

{u' —u—ie) k 

1 



fc-1 (u'-u-ie) k - 1 
k-l 



In 2U(w' -u-ie) 
In 2H{u — u — ie) 



+ In 
+ ln 



-2H(u' - u-ie) } 



(4.13) 



—2H(u — u — ie) 

k-l 



+ 



k-l 



i(f) " [H*He) + ^] - ^(jJ^) ~ {^[2(1-^)]+^+^} , 
where $(z, s,v) is the special function [14]: 



00 



n=0 



(v + nY 



(4.14) 



and the arrow indicates the leading infinitesimal e contributions. When no logarithms are 
present the integrals are simple enough that we give only the small e forms: 



(i) 



H- 



du' 



w — u — ie 



- ln(l - Hu) - \n{He) + £tt , 
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(4.15) 



(iv) 
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(4.16) 
(4.17) 
(4.18) 



{u! — u — ie) k 

We can now complete the reduction of the contribution from the i = 1, j = 10 vertex 



pair: 



k a H 2 



1,10 " 26^4 
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-61n 2 (ifw) + 



^ - 12 In 2 



ln(ffu) + 19(1 - Hu) + § 7T 



{( 

+ 12 [C(2) - Hu$(Hu, 2, 1)] - ^ In 2(1 - ifu) 
+ f ln(2tf e ) + 61n 2 (2tf e )}-L + ^(t^)^ 

.(l=Hu) + ?]^ + I["(t 



3^2 



-6 In 2 



2(1 - ifu) 



15 
~8~ 



(4.19) 



It is worth noting that the coefficient of the leading term for late times 6 u 4 In 2 (.Hit) 

- is opposite to that of the double logarithmic ultraviolet divergence - +6 ^ 4 ln 2 (2ife). 
This is in sharp contrast to 3-3-3 contributions such as (3.47a) where an undifferentiated 
propagator logarithm is not involved. Then the leading infrared contribution has the same 
sign as the double logarithmic ultraviolet divergence. 

The same sign phenomenon of (3.47a) was explained in sub-section 2.7. The parameter 
e has dimensions of length so it can only appear in the combinations He or eu . The 
dimensionality of the conformal time integrands shows that they converge even if the 
upper limit, -ff -1 , diverges, so the only factors of ln(H) can come from the single possible 
undifferentiated propagator logarithm. When no such logarithm is present any double 
logarithmic ultraviolet divergence must take the form: 

\n 2 (eu- 1 ) = ln 2 (He) - 2 ln(He) ln(Hu) + ln 2 (Hu) . (4.20) 

One can understand why the opposite relation prevails in all 4-3 contributions by the 
need to avoid overlapping divergences. The 4-3 diagram contains only two vertices so its 
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ultraviolet divergences must be primitive; there cannot be any non-local sub-divergences 
such as \n(He) ln(Hu). The possibility of In 2 (He) is ruled out because no factors of 
ln(H) can come from the limits of integration; the only ln(H) can come from the single 
possible undifferentiated propagator logarithm. So we are left with just one dimensionally 
consistent double ultraviolet logarithm: 

ln 2 (#e) - ln 2 (Hu) = ln 2 (e) + 2 ln(e) ln(H) - 2 \n(H) ln(u) - ln 2 (u) . (4.21) 

Although the tensor algebra and derivatives of the 4-3 diagram are only marginally 
simpler than that of the most complicated 3-3-3 diagram (2a), the number of terms they 
produce is much smaller. There are just 39 integrands which contain ln(H 2 w 2 ) and only 
57 which have no logarithm. We have also seen that the integrals are simple enough for 
exact results to be obtained. Although we did this, we report just the leading results here: 

, , n A H 2 f 13 ln 2 (M „/\n.(Hu)\\ , * nn s 

, , k A H 2 f ln 2 (M „An(Hu)\) 

<*»(«) = { +2 4^ + (^)} • (4226) 



5. Epilogue 

We can now assemble the results from the various diagrams of Fig. 2. Combining (3.51) 
with (4.22) gives the following coefficient functions for the amputated 1-point function 
(2.19): 

a(u) = H 2 (^—J |( _+_+___) \n 2 (Hu)+0(\n(Hu)) j+0(« 6 ) (5.1a) 

<") = H ~ 2 (S) 4 {( +1 T ! -T" 112 + 8 ) \AHu) + 0(HHu)) }+0(„ 6 ) (5.16) 
The four numbers in each of these expressions represent the respective contributions from 
diagrams (2a), (2b), (2c) and (2d). As explained in sub-section 2.9, diagram (2e) can only 
contribute a single factor of ln(Hu) and diagram (2f) cannot contribute any. 
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Attaching the retarded Green's functions discussed in sub-section 2.5 gives the coeffi- 
cient functions for the full 1-point function (2.20): 

2J -j In 3 (Hu)+0 (in 2 (Huj) |+0(k 6 ) (5.2a) 
52 + ll) ln 2 (Hu) + o(ln(Hu)^ \+0(k 6 ) . (5.26) 



^-(s) 4 {(- 
C «-(S) 4 {C 



7180 2416 1280 208\ 
f— - ^ + 
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81 



27 



319 49 



None of the four coefficient functions given above is free of gauge dependence. The physical, 
gauge independent observable is the combination (2.24) which gives the effective Hubble 
constant. Our result for it is: 

4309 1649 172 



+ 
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9 



+ -) (Htf + 0(Ht) 



Hll- 



+ O(k ) 

eH) 4 [> 2 +°H 



(5.3a) 
(5.36) 



Although we have given the contribution from each diagram in (5.3a), only the sum (5.3b) 
is really physical and gauge independent. Still, it is worth pointing out that the pure 
graviton diagrams — (2a) and (2d) — act to slow inflation while the ghost diagrams - 
(2b) and (2c) — contribute in the opposite sense. 

The physical significance of our result has been discussed elsewhere [4] but two points 
should be mentioned here. First, we get a reduction in the effective Hubble constant due 
to the negative energy of the gravitational interaction between the zero point motions of 
gravitons.* This tendency ought to persist to all orders in perturbation theory. Second, 
recall from sub-section 2.8 that at most I infrared logarithms can appear at I loops. This 
means that the £-loop contribution to the bracketed term in (5.3b) is at most: 



-# (kH) 2£ (Ht) 



(5.4) 



* The induced stress-energy is that of negative vacuum energy, to leading order, because causality restricts the 
gravitational interaction to the constant Hubble volume. The induced energy density is therefore independent 
of the true volume of the inflating universe, so the total energy is just E = Vp and the pressure is p = 
-dE/dV = -p. 
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The bound is quite likely to be saturated because one can form higher loop graphs by 
attaching the 2-loop tadpole we have computed. But one in any case obtains the following 
estimate for the number of e-foldings needed to end rapid inflation: 

N = Ht~ (kH)- 2 . (5.5) 

This is > 10 12 for inflation at the GUT scale or below, which is more than enough to 
explain the homogeneity and isotropy of the observed universe. If the unknown numerical 
coefficients in (5.4) fall off less rapidly than then the end of rapid inflation is likely 
to be quite abrupt because all orders will become strong at once. This augurs well for 
reheating. 

We turn now to a discussion of accuracy. The strongest checks on the consistency of our 
basic formalism and the accuracy of our implementation are provided by the one-loop self- 
energy. This quantity was relatively simple to obtain [15] because technical considerations 
compelled us to compute the inner loops of diagrams (2a) and (2b) before doing the 
outer loop contractions and acting the outer loop derivatives. We checked that the one- 
loop self-energy has the appropriate reflection symmetry, that it obeys the Ward identity, 
and that it agrees, for H — > at fixed t, with the flat space self-energy obtained by 
Capper [16] in the same gauge. Since we used only partially symmetrized vertices and 
then interchanged lines where necessary, reflection symmetry is a non-trivial test of our 
programs for the tensor algebra. It also checks the programs for taking derivatives because 
the order of differentiation breaks manifest reflection symmetry. The Ward identity tests 
the apparatus of gauge fixing, our solutions for the ghost and graviton propagators, and 
the 3-point interaction vertices. It also provides a powerful independent check of the tensor 
algebra and derivative programs. In addition to giving further independent tests for all 
these things, the flat space limit checks the overall factor and the sign. 

One of the major complications in computing diagram (2a) was that the number of 
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intermediate expressions becomes prohibitive if one attempts to perform the tensor algebra 
all at once, even for a single "trio" of vertex operators.* This is why we first did the 
inner loop tensor algebra and acted the inner loop derivatives for each pair of inner loop 
vertices, then summed the results and projected the lengthy total onto a basis formed from 
products of the 79 independent 4-index objects and the 10 possible contracted outer leg 
derivatives. Each basis element with a non-zero coefficient was contracted into the outer 
loop propagators and the outer vertex operators, and the various outer loop derivatives 
were acted. Then the inner loop coefficients were multiplied and the results summed. 

The cumbersome nature of this procedure caused us much anxiety and we devised 
a number of ways to check it. First, a program was written to check the inner loop 
projection by simply summing the product of each basis element with its coefficient and 
then comparing with the original expression. Second, the outer loop tensor algebra and 
derivative programs were based on the same scheme as those which were so effectively 
checked by the Ward identity, and of course we used the same stored expressions for the 
vertices and propagators. Third, it was not too time consuming to perform all the tensor 
algebra at once for a single "trio" of vertex operators. This was done for the case where 
the vertex operator at x^ 1 is #10 in Table 3, and where those at x and x"^ are both 
#41. The derivatives were then acted (using a different program from the usual one) and 
the result compared with what our standard programs give. The same calculation was 
performed by hand as an additional check. 

Diagram (2b) was computed using the same programs as for (2a), so its accuracy is 
checked by that of (2a) . An additional test was provided by the fact that no undifferenti- 
ated propagator logarithms can come from ghost lines.* We checked that the integrands 

The problem is that expanding the four propagators gives 6 4 = 1296 terms, each of which involves a 
contraction over 16 indices. 

To see this note that the 10 vertices of Table 2 have either a a factor of <9 2 or a factor of t a ' 2 — which 
accesses only the normal part of the propagator. 
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for (2b) are indeed free of ln(if 2 y 2 ), and that those of diagram (2c) are free of either 
\n(H 2 w 2 ) or ln(#V). 

Our greatest worry was the complicated integrations of 3-3-3 diagrams. We checked 
these a number of ways. First, we did a large number of examples by hand. Second, 
the two authors wrote independent integration schemes. The first of these did each radial 
integration independently by the method of contours and changed variables as necessary so 
that only the final conformal time integration needed to be expanded. This program also 
computed the four "±" variations separately. The second program did the radial integrals 
by differentiating generating functions. The conformal time integrations were carried out 
in the standard order, expanding whenever necessary. This second program also added the 
"±" variations of the integrands before evaluating the integral. We ran the two programs 
for each of the thousands of terms and then scanned the results to make sure they agreed. 

One of the strongest checks on the accuracy of the 3-3-3 diagrams derives from our 
proof in sub-section 2.8 that each "trio" of vertices for each of the 3-3-3 diagrams must be 
free of triple log terms — tt -4 ln 3 (ifw). As mentioned there, and as we saw explicitly in 
(3.49), this is not true for each of the integrands which emerge from step 3 of our reduction 
procedure. Yet when the results from integrating the thousands of distinct integrands were 
summed, the triple logs cancelled for each diagram. We also checked that our programs 
show this cancellation separately for the previously mentioned 10-41-41 vertex triad. 

We fretted long about how to check the 4-3 diagram (2d). Of course it is partially 
checked by the 3-3-3 diagrams by virtue of the fact that all use the same stored expression 
for the graviton propagator, and since their tensor algebra and derivative action programs 
are based on the same scheme. One of our big worries was the different vertices. We checked 
the fully symmetrized 3-point vertex by having the computer symmetrize the partially 
symmetrized vertex of Table 3. We also wrote a program in which the computer produces 
fully symmetrized 3-point and 4-point vertex operators from the respective interaction 
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Lagrangians (3.6) and (4.1). These were compared directly with our stored expression for 
the symmetrized 3-point vertex; of course we had to (computer) symmetrize our 4-point 
vertex before comparing it. 

As an additional check on the 4-3 diagram, we compared the output from our programs 
with a hand computation of the result for vertex operator #1 at x^ and #10 at x 1 ^ 1 . It is 
also worth pointing out that the 4-3 integrations are vastly simpler than those of the 3-3-3 
diagrams. In fact only the single term (4.11) can contribute at leading order. 

Finally, there is the fact that the result is plausible. Quantum gravity is not on-shell 
finite at one loop when the cosmological constant is non-zero [17], so one has to expect 
ultraviolet divergences of the form ln 2 (e) at two loops. Since e has dimensions of length, 
logarithms of it must come in the form ln(ett _1 ) and ln(He). We showed in sub-section 2.7 
that only a single factor of ln(if) can occur, so there must be at least one hi(u), and we can 
think of no reason why two should not occur. The reason why these infrared logarithms 
act to slow inflation is that they represent the negative gravitational interaction energy 
between the zero point motions of gravitons. At the risk of putting too much faith in 
gauge-dependent results one can even understand that pure graviton diagrams should act 
to slow inflation, and that the ghost diagrams should diminish this effect.* 

There is also a good physical interpretation for infrared logarithms. The up to two 
powers of ln(Hu) that derive from integrating factors of it and u" represent the 
invariant volume of the past lightcone as viewed from the observation point [2]. The 
single ln(Hu) which can come from an undifferentiated propagator logarithm represents 
the ever-increasing correlation of the free graviton vacuum in an inflating universe. It has 
long been known that the assumption of correlated de Sitter vacuum over an infinite surface 
of simultaneity leads to a divergence in the propagator. (This was first proved by Allen 

The near cancellation between the two classes is perhaps explicable from the fact that the gauge fixed 
graviton field carries eight unphysical modes and only two physical ones. At two loops this means roughly 
10 2 = 100 mode pairs of which the ghost loops must remove all but 2 2 = 4. 
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and Folacci [18]. See also [7] and the references cited therein.) Infrared logarithms are the 
causal manifestation of this effect when one approaches an infinite surface of simultaneity 
by inflating a finite patch of correlated de Sitter vacuum. 

In sum, the procedure we used should work, it gives a reasonable result, and every care 
has been taken to ensure accuracy. However, we do not wish this discussion to convey a false 
sense of infallibility. We have been schooled in humility by the disconcerting experience of 
detecting errors even after the final one seemed to have been expunged. It must be stressed 
that this was an enormously complicated piece of work. Only one other two-loop result 
has ever been obtained in quantum gravity [19], and it was confined to the ultraviolet 
divergent part of the on-shell effective action for zero cosmological constant. We feel very 
strongly that computer calculations on this scale should be regarded as experiments whose 
results require independent verification before they can be completely trusted. We would 
be happy to cooperate with anyone wishing to undertake even a partial check. 

We conclude by briefly discussing a perturbative issue which demands further study: 
the case of a negative cosmological constant. This is interesting in its own right and because 
it may have relevance to the period after rapid inflation has ended, when an energetically 
favored phase transition would be expected to generate a negative effective cosmological 
constant. Two qualitative questions are of great importance: are there strong infrared 
effects from quantum gravity? and, do they tend to resist the contraction of spacetime? 

Careful consideration of the effect for positive cosmological constant leads to the con- 
clusion that it derives from the combination of three features: 

(1) Propagators which do not oscillate or fall off over large temporal separations; 

(2) An interaction of dimension three; 

(3) The fact that the invariant volume of the past lightcone increases without bound. 

The last two are certainly true as well for a negative cosmological constant. In fact the 
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causal structure of anti-de Sitter space allows one to access spatial infinity after only a 
finite amount of time. The non-trivial issue is the propagator. The usual prescription for 
defining what happens at spatial infinity is to impose reflective boundary conditions [20]. 
These are enforced mathematically by a negative image source at the antipodal point, so 
they make propagators fall off too rapidly to give a big infrared effect. 

We believe that reflective boundary conditions are reasonable for Euclidean anti-de 
Sitter space because its antipodal points are not part of the manifold. However, the an- 
tipodal points of the Minkowski signature formalism are on the manifold. In this case the 
use of reflective boundary conditions implies the physical absurdity that every source of 
stress energy has an antipodal anti-source. It is not reasonable to suppose a man is for- 
bidden to shake his fist (which generates gravitational radiation) without the cooperation 
of someone of the other side of the Universe. 

The more sensible boundary conditions seem to be transmissive, in which information 
is allowed to flow to spatial infinity without hindrance. These do not lead to a globally 
well-posed initial value problem, but they do make physical sense locally. The propagator 
associated with this condition does not fall off rapidly enough to preclude a large infrared 
effect. It should also be noted that arguments in the literature about the stability of super- 
gravity or superstrings on an anti de Sitter background are based on reflective boundary 
conditions.* 

It is very much more difficult to determine whether quantum gravity makes a negative 
A universe collapse slower or faster. Local considerations compel us to the view that 
gravitational interaction energy should still be negative. However, the fact that the sign 
of the dimension three coupling changes makes it hard to say what this does. (Note that 
the 1-point function is odd in the 3-point vertex.) An additional complication is that we 

can no longer count on the induced stress tensor to be that of pure vacuum energy. That 
We are indebted to I. Antoniadis for bringing this point to our attention. 
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it came out this way for the case of positive cosmological constant derives from the finite 
causal horizon of de Sitter space, which means that the gravitational interaction energy 
density must be independent of the true volume. This is not true for anti-de Sitter space. 
And it is worth recalling that inflation is driven by the negative pressure of a positive 
cosmological constant; the positive energy density serves as a drag on expansion. So it is 
conceivable that negative A quantum gravity induces a negative energy density, while still 
resisting contraction by virtue of generating less than an equal and opposite pressure. 



ACKNOWLEDGEMENTS 

We wish to give special thanks to T. J. M. Zouros for the use of his SUN SPARC station, 
and to S. Deser for his encouragement and support during this project. We have profited 
from conversations with I. Antoniadis, C. Bachas and J. Iliopoulos. One of us (RPW) 
thanks the University of Crete and the Theory Group of FO.R.T.H. for their hospitality 
during the execution of this project. This work was partially supported by DOE contract 
86-ER40272, by NSF grant 94092715 and by EEC grant 933582. 



63 



REFERENCES 

[1] N. C. Tsamis and R. P. Woodard, Phys. Lett. B1993 (1993) 351. 

[2] N. C. Tsamis and R. R Woodard, Ann. Phys. 238 (1995) 1. 

[3] S. Deser and L. F. Abbott, Nucl. Phys. B195 (1982) 76. 
P. Ginsparg and M. J. Perry, Nucl. Phys. B222 (1983) 245. 

[4] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B474 (1996) 235. 

[5] N. C. Tsamis and R. P. Woodard, Commun. Math. Phys. 162 (1994) 217. 

[6] N. C. Tsamis and R. P. Woodard, Phys. Lett. B269 (1992) 269. 

[7] N. C. Tsamis and R. P. Woodard, Class. Quantum Grav. 11 (1994) 2969 

[8] J. Schwinger, J. Math. Phys. 2 (1961) 407; Particles, Sources and Fields (Addison- 
Wesley, Reading, MA, 1970). 

[9] L. H. Ford, Phys. Rev. D31 (1985) 710. 

[10] A. D. Dolgov, M. B. Einhorn and V. I. Zakharov, Phys. Rev. D52 (1995) 717. 

[11] B. S. DeWitt, Phys. Rev. 162 (1967) 1239. 

F. A. Berends and R. Gastmans, Nucl. Phys. B88 (1975) 99. 

[12] S. Wolfram, Mathematica, 2nd edition (Addison- Wesley, Redwood City, CA, 1991). 

[13] R. Mertig, Guide to FeynCalc 1.0, University of Wiirzburg preprint, March 1992. 

[14] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, 4th edition 
(Academic Press, New York, 1965) pp. 1075 - 1076. 

[15] N. C. Tsamis and R. P. Woodard, Phys. Rev. D15 (1996) 2621. 

[16] D. M. Capper, J. Phys. A13 (1980) 199. 

[17] S. Deser and P. van Nieuwenhuizen, Phys. Rev. D10 (1974) 401. 

64 



[18] B. Allen and A. Folacci, J. Math, Phys. 32 (1991) 2828. 

[19] M. Goroff and A. Sagnotti, Phys. Lett. B160 (1986) 81; Nucl. Phys. B266 (1986) 
709. 

[20] S. J. Avis, C. J. Isham and D. Storey, Phys. Rev. D18 (1978) 3565. 



65 



